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INTRODUCTION 

BACKGROUND 

DMM 

A  particular  dynamic  design-anaiysis  method.  DDAM,  has  been 
applied  to  selected  equipment  on  Navy  ships  since  I960.  The  DDAM  is 
based  on  measurements  of  the  responses  of  shipboard  equipment  to 
controlled  tests  against  ships  with  underwater  explosions. 
Measurements  from  a  variety  of  ships  showed  that  the  responses  of 
equipment  to  the  shock  from  an  underwater-explosion  attack  could 
be  approximated  fairly  well  by  considering  only  the  location  of  the 
equipment  on  the  ship,  the  weight  of  the  equipment,  and  its  natural 
frequency  of  vibration.  Collected  measurements  were  fitted  to  simple 
formulas  depending  on  location,  weight,  and  frequency  that  could  be 
used  to  check  the  ability  of  equipment  to  resist  shocks  from 
underwater  explosions. 

The  DDAM  can  be  looked  at  in  modern  terms  as  an  approximate 
method  of  conducting  a  substructure  analysis  of  a  complex  dynamic 
system.  In  the  DDAM  all  of  the  effort  in  the  analvsis  is  directed 
toward  one  particular  substructure,  representing  the  equipment  of 
interest.  The  main  structure,  consisting  of  ail  the  rest  of  the  ship,  is 
taken  as  a  standard  structure  whose  properties  do  not  need  to  be 
specially  determined.  The  experimental  measurements  of  the  shock 
that  the  ship  can  transmit  to  equipment  of  different  weights  and 
frequencies  during  an  underwater-explosion  attack  take  the  place  of 
an  analysis  of  the  ship  s  structure 

The  effort  of  the  analvst  in  applying  the  DDAM  is  limited  to 
determining  the  frequencies  and  weights  associated  with  one 
particular  item  of  equipment.  The  formulas  from  the  DDAM  then 
provide  the  analyst  immediately-  with  the  responses  that  would  be 
expected  at  each  weight  and  frequency  if  the  equipment  were 
mounted  at  a  particular  location  on  a  standard  ship  during  an  attack 
against  the  ship  with  an  underwater  explosion  The  analvst  performs. 
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in  effect,  a  substructure  analysis  without  having  to  analyze  the  mam 
structure  or  consider  location  of  the  equipment  bevond  the  broad 
categories  given  in  the  formulas 


IfiDAM 

Calculations  for  the  DDAM  are  based  on  the  assumption  that  the 
shock  from  an  underwater  explosion  produces  a  translation  of  the  ship 
in  a  single  direction,  with  all  of  the  supports  for  the  equip ment 
moving  identically.  This  is  generally  a  good  approximation  for 
compact  equipment  or  for  explosions  that  take  place  in  the  water 
some  distance  from  the  ship 

An  extension  of  the  DDAM,  referred  to  as  IDDAM  (Impact 
DDAM)  was  developed  in  1977  in  response  to  a  request  for  a  dynamic 
design-analysis  method  that  could  be  applied  to  the  case  of  the  impact 
of  exercise  (nonexploding)  torpedoes  against  the  hull  of  a  submarine 
Here  the  shock  applied  to  the  several  supports  for  a  particular  item  of 
equipment  attached  to  the  hull  could  vary  greatly,  from  a  severe 
shock  to  a  support  directly  in  way  of  the  impact  to  shocks  of  smaller 
and  smaller  severity  with  increase  in  distance  from  the  point  of 
impact. 

Determining  the  response  of  a  structure  to  independent  motions 
of  multiple  supports  involved  only  a  minor  elaboration  of  DDAM 
However,  the  weight  of  the  structure  had  to  be  distributed  among  the 
separate  supports  in  a  self-consistent  fashion  in  order  to  allow  the 
simple  formulas  involving  weight  and  frequency  to  be  applied  When 
all  of  the  supports  move  identically,  as  in  DDAM.  the  full  weight  reacts 
against  the  common  motion,  but  when  the  motions  are  different  the 
weight  appears  as  a  combination  of  direct  and  cross  terms  among  the 
supports,  adding  appreciably  to  the  complexity  of  the  method. 

REVISION  OF  THE  IDDAM 


The  present  report  describes  a  new  version  of  the  IDDAM 
referred  to  here  as  RIDDAM  'Revised  IDDAM).  The  RIDDAM  simplifies 
the  method  ol  calculating  the  effective  weight  of  a  structure  with 
multiple  supports 
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The  original  DDAM  was  developed  at  a  time  when  structural  - 
analysis  computations  were  regularly  done  using  pencil,  paper,  and  a 
desk  calculator.  Its  procedures  have  been  carried  over  to  the  present 
IDDAM. 

First,  turn  on  your  computer.  In  the  following,  the  equations  for 
the  IDDAM  have  been  written  using  the  matrix  notation  that  is 
common  for  computer-aided  analysis,  on  the  assumption  that  a 
computer  rather  than  a  desk  calculator  will  be  used  to  analyze  the 
structure.  Appendix  A  gives  some  of  the  elements  of  matrix  algebra 
for  readers  who  are  not  entirely  familiar  with  matrix  notation 
Appendix  B  lists  a  computer  program  that  can  be  used  to  analvze 
simple  structures  on  any  computer  that  has  a  compiler  or  interpreter 
for  the  BASIC  computer  language. 

The  rewritten  equations  for  the  IDDAM  produce  a  more  compact 
notation  and  simpler  calculations  but  make  no  change  in  its 
procedures. 


The  RIDDAM  does  make  a  substantial  change  in  procedures  bv 
the  method  it  uses  to  determine  effective  weight,  however 

In  both  the  DDAM  and  the  IDDAM  the  weight  of  a  normal  mode 
of  vibration  is  considered  as  reacting  against  supports  only  to  the 
extent  that  the  supports  are  driven  by  a  shock  input.  The  weight  of  a 
mode  thus  varies  with  direction  of  the  input  in  the  DDAM  and  varies 
with  the  distribution  of  shock  severities  among  the  supports  in  the 
IDDAM. 

The  method  of  determining  the  weight  of  a  mode  in  the 
RIDDAM  presumes  that  the  mode  reacts  against  all  of  the  supports 
whenever  it  responds  to  a  shock  delivered  through  any  of  them  The 
effective  modal  weight  in  the  RIDDAM  is  thus  independent  of  the 
particular  shock  that  is  applied  Calculations  are  not  only  simplified 
but  reduced  to  the  more-rational  basis  of  assuming  that  the  structure 
reacts  against  all  of  its  supports  whenever  it  is  forced  to  undergo 
accelerations  from  anv  cause 


PROCEDURE 


The  equipment  (substructure)  is  modeled  as  a  linear,  elastic, 
and  undamped  structure.  Details  of  the  model  depend  on 
requirements  of  the  computer  program  that  is  used. 


Frequencies  are  found  by  doing  a  standard  normal -mode 
analysis  of  the  model  with  fixed  interfaces  at  the  points  of  attachment 
to  the  ship. 


The  normal-mode  analysis  provides  a  set  of  factors  showing 
how  the  motion  of  each  support  participates  in  producing  a  response 
of  each  of  the  normal  modes 


The  critical  feature  of  both  the  DDAM  and  the  IDDA.V1  is  the 
assignment  of  an  effective  mass  to  each  normal  mode  of  vibration 
The  effective  mass  is  used  as  a  measure  of  the  reactive  force  that  the 
mode  applies  to  the  supports  for  the  equipment  when  the  equipment 
is  accelerating  in  response  to  a  shock.  In  the  R I  DDAM  this  effective 
mass  is  obtained  simply  by  adding  mass-normalized  participation 
factors  by  absolute  value  and  squaring  the  sum. 


Input? 


As  for  the  DDAM,  the  inputs  for  the  IDDAM  are  formulas 
showing  the  response  of  a  normal  mode  of  vibration  as  a  function  of 
its  weight,  frequency,  and  location  on  board  ship  In  the  DDAM  all  of 
the  supports  are  assumed  to  move  together  in  one  particular  direction 
The  IDDAM  relaxes  this  assumption  by  providing  different  inputs  for 
different  supports.  The  RIDDAM  uses  the  same  inputs  as  the  IDDAM 
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Responses  of  the  modes  are  combined  to  estimate  forces, 
stresses,  strains,  and  deformations  throughout  the  model  of  the 
equipment.  The  stresses  and  deformations  at  critical  locations  are 
compared  with  allowable  values  to  assess  the  ability  of  the  equipment 
to  withstand  a  particular  shock.  Deformations  produced  by 
differences  in  the  motions  of  redundant  supports  are  included  in  both 
the  IDDAM  and  the  RIDDAM. 

OUTLINE 

The  report  begins  with  a  brief  overview  of  the  theoretical  basis 
for  the  normal-mode  method  used  in  the  RIDDAM.  It  continues  with  a 
discussion  of  some  of  the  practical  compromises  that  are  necessary  to 
keep  its  implementation  simple,  and  a  description  of  some  of  the 
problems  that  the  compromises  generate.  There  are  some  worked-out 
examples  of  analyses  of  simple  structures  by  the  RIDDAM.  The  object 
is  to  supply  enough  information  so  that  an  analyst  can  apply  the 
RIDDAM  to  his  own  item  of  equipment  by  following  the  discussions 
and  examples  given  here. 


THEORY 


MATHEMATICAL  MODEL 


Most  computer  programs  for  dynamic  analvsis  represent  the 
structure  by  a  matrii  equation  of  the  form 

Ma  ♦  K  x  -  Fu.  <  I ) 

where  alt)  is  a  time-varying  column  vector  of  accelerations  at  selected 
points  (nodes)  on  the  structure  and  x(t)  is  a  column  vector  of 
displacements.  The  symmetric  matrices  M  and  K  represent  the  mass 
and  stiffness  of  the  structure.  The  column  vector  u(t)  consists  of 
displacements  of  the  supports  for  the  structure;  the  rectangular  matrix 
of  support  stiffnesses  F  converts  displacements  of  the  supports  into 
forces  applied  to  the  coordinates  of  the  structure  The  equation 
balances  these  forces  against  the  forces  Ma  and  Kx  from  the  mass  and 
stiffness  of  the  structure. 

The  motions  represented  by  Equation  1  can  be  translations  in 
any  direction  or  rotations  about  any  axis,  with  elements  of  M.  K.  and  F 
being  mass,  mass  moment  of  inertia,  force  per  deflection,  or  moment 
per  angle.  It  is  assumed  that  all  the  translations  and  rotations  are 
small,  the  stiffnesses  are  constant,  and  there  is  no  damping  in  the 
structure. 

Discussions  here  are  framed  in  terms  of  translations 
(displacements  and  forces)  with  the  understanding  that  rotational 
motions  (angles  and  moments)  are  also  included 


Some  methods  of  structural  analysis  are  based  on  other  than 
the  force-balance  equation.  For  example,  the  displacement-balance 
equation 


K-l  M 
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has  been  popular  for  use  with  a  desk  calculator.  Here  K1  'the  inverse 
of  the  stiffness  matrix)  is  a  matrix  of  influence  coefficients  that  show? 
the  displacements  of  all  the  points  on  the  structure  when  unit  force  is 
applied  to  one  of  the  points.  The  support  influence  matrix  K~'F 
represents  the  position  of  the  structure  produced  by  the  motion  of  its 
supports.  The  equation  balances  these  equilibrium  displacements 
against  the  dynamic  displacements  that  include  forces  from  the 
accelerations  of  the  masses. 

The  displacement-balance  equation  is  not  well  suited  to 
computer-aided  analyses  because  the  matrix  K*M  is  full  of  nonzero 
elements  and  is  not  symmetric.  A  computer  needs  less  memorv  and 
can  run  faster  using  the  matrices  M  and  K.  both  of  which  are  sparse 
and  symmetric. 

Discussions  here  are  in  terms  of  the  force-balance  equation  as 
the  equation  most  common  in  computer  programs  for  dynamic 
structural  analysis 

NORMAL  MODES  OF  VIBRATION 


Many  computer  programs  begin  a  normal- mode  analysis  by- 
generating  a  square  or  rectangular  matrix  Q  that  makes 

0T  M  Q  -  I  1 3 1 

be  a  unit  matrix.  Normalizing  on  the  mass  matrix  in  this  way- 
produces  a  simple  and  efficient  analysis,  and  has  been  assumed  in  the 
following  discussion.  A  later  paragraph  will  consider  the  alterations 
that  are  necessary  when  other  methods  of  normalization  are  used 

Diagonalizing  the  stiffness  matrix 

The  matrix  Q  can  also  be  chosen  to  make 


QTKQ  =  p2 
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be  a  diagonal  matrix.  There  are  many  different  numerical  methods 
for  generating  a  Q  that  will  satisfy  Equations  3  and  4  simultaneously. 
Which  is  best  depends  on  the  size  of  the  problem,  the  distribution  of 
the  nonzero  elements  in  M  and  K,  and  how  many  elements  of  Q  are  to 
be  calculated.  Manv  programs  offer  a  choice  of  methods  and  may  also 
have  ways  of  decreasing  the  sizes  of  the  structural  matrices  to  speed 
up  the  calculation. 

Normal-mode  equations 

If  the  matrix  Q  is  used  to  transform  the  accelerations  and 
displacements  by 

a-Qb.  :  >  Q  y,  i  3 ' 

and  if  Equation  1  is  premultiphed  by  QT  there  results 

b  ♦  p2  y  _  qT  p  u  1 6 ' 

Because  P2  is  diagonal,  each  element  of  the  acceleration  bit)  and 
displacement  y(tl  on  the  left  side  represents  an  independent  oscillator 
with  a  fixed-base  angular  frequency  given  by  the  square  root  of  its 
element  of  P2  The  responses  of  the  oscillators  can  be  converted  back 
to  motions  of  the  structure  bv  Equation  5 

Standard  form 

Equation  6  can  be  reduced  to  standard  form  for  each  oscillator 
by  writing  it  as 

b  »  P2  y  »  p2  RT  u  [7  i 

where 

RT  .  p  2  qT  F  1 8  I 


is  a  rectangular  matrix  called  the  participation  factor  The  matrix  RT 
converts  the  displacements  uu1  of  the  supports  to  a  generalized 


displacement  representing  the  motion  of  the  base  of  each  one  of  the 
normal-mode  oscillators.  The  participation-factor  matrix  has  as  many 
columns  as  there  are  separate  supports  for  the  structure,  and  as  many 
rows  as  there  are  normal  modes. 

Each  of  the  normal  modes  in  Equation  7  responds  as  if  it  were  a 
single  oscillator  supported  from  a  base  whose  displacement  is  given 
by  one  of  the  elements  of  RTu(t).  That  is,  if  y(t)  is  the  displacement 
and  b(t)  is  the  acceleration  of  one  of  the  normal  modes,  the  equation 
for  its  response  is 

b  +  p2  y  =  p2  fT  u ,  ( 9 ) 

where  rT  is  the  row  from  RT  that  corresponds  to  the  elements  y»i)  and 
bit),  and  P2  is  the  square  of  the  angular  frequency  of  the  mode. 

MODAL  MASS 

What  the  modal  mass  does 

A  force  is  required  to  produce  an  acceleration  of  a  real  structure 
that  has  mass.  This  force  reacts  against  the  supports  for  the  structure 
and,  in  all  practical  cases  where  the  supports  are  attached  to  some 
other  structure,  produces  motions  of  the  supporting  structure. 

The  displacements  uiti  of  the  supports  need  to  be  considered  as 
consisting  of  two  components  One  component  is  the  motion  »or  lack  of 
motion)  that  would  be  expected  if  the  structure  were  not  in  place  The 
other  component  is  the  motion  produced  bv  the  reaction  force  as  the 
structure  accelerates. 

If  the  structure  is  analyzed  into  its  normal  modes,  each  normal 
mode  needs  to  be  assigned  a  mass  so  its  reaction  force  can  be 
determined. 

A  simple  example  of  modal  mass 

Suppose  a  single  mass  M  is  supported  by  a  set  of  parallel 
springs  having  a  total  stiffness  k  The  equation  of  motion  can  be 
written 


M  a  *  k  i  -  fT  u, 


where  fT  is  a  row  vector  whose  elements  are  the  stiffnesses  of  the 
individual  springs  and  u(t)  is  a  column  vector  of  the  displacements  at 
the  points  where  each  spring  is  attached  to  some  other  structure 

Follow  the  method  of  analysis  given  in  the  preceding  sections  as 
if  the  mass  M  and  the  stiffness  K  were  matrices  with  just  one  row  and 
one  column.  First  define  a  modeshape  Q  that  makes 

QMQ  •  1.  '  1  1 1 

Notice  that  Q  is  the  inverse  square  root  of  M.  and  that 

QKQ  -  K/M  -  P2  '  12  i 

is  the  square  of  the  natural  frequencv  of  the  mass  as  supported  by  all 
the  springs. 

Continue  to  find  the  participation  factor 
rT  -  p-2  q  fT  1 1 3 1 

and  reduce  it  to 

rT  *  (K/MH  M-l/2  fT  -  vM  fT/k.  (Mi 

Summing  the  elements  of  the  participation-factor  vector  gives  the 
square  root  of  M.  and  squaring  the  sum  recovers  M  as  the  mass  of  the 
single  mode  of  the  structure 

The  same  procedure  can  be  applied  to  Equation  9  to  determine 
an  effective  mass  for  the  mode  whose  acceleration  is  b(tl  That  is  add 
the  elements  of  the  row  from  the  participation-factor  matrix  and 
square  the  sum  to  find  the  mass  that  must  be  accelerated  when  the 
structure  responds  with  motions  in  one  particular  normal  mode  of 
vibration. 


The  advance  of  the  DDAM  was  in  pointing  out  that  the  effect  of 
reaction  forces  on  a  particular  structure  can  be  found  experimentally 
by  mounting  substructures  to  it  and  measuring  the  responses  of 
substructures  having  different  masses  and  frequencies.  The  results 
from  the  tests  can  then  be  applied  to  predict  the  responses  of  modes 
of  other  items  mounted  to  the  structure  provided  their  masses  and 
frequencies  are  known.  The  supporting  structure  thus  does  not  need 
to  be  analyzed  in  order  to  estimate  the  shock  response  of  items 
attached  to  it. 

Some  approximation  is  inherent  in  this  assumption,  especially  if 
responses  measured  for  one  particular  structure  are  applied  to 
different  but  similar  structures.  Supporting  structures  can  differ 
because  of  other  substructures  mounted  near  the  substructure  of 
interest,  or  even  as  a  result  of  other  modes  of  the  substructure  being 
analyzed.  However,  the  advantage  of  being  able  to  estimate  the 
response  of  a  particular  item  of  equipment  without  having  to  make  a 
detailed  analysis  of  the  supporting  ship's  structure  makes  the  method 
attractive  despite  the  potential  errors. 


Squaring  the  sum  of  the  elements  in  a  row  of  the  influence- 
coefficient  matrix,  as  shown  in  the  simple  example  above  will  give  an 
effective  mass  for  the  mode  if  all  the  reactions  are  in  the  same 
direction,  as  they  were  for  the  parallel  springs  in  the  example  If 
some  of  the  reactions  are  positive  and  some  are  negative,  the  squared 
sum  will  give  an  effective  mass  producing  a  net  force  on  all  of  the 
supports  taken  together.  This  is  the  effective  mass  used  in  the  DDAM. 
where  differences  in  the  motions  of  different  supports  are  not 
considered. 

If  each  point  of  support  is  taken  as  attached  to  a  separate 
structure,  the  reactions  will  produce  independent  motions  of  each 
support  in  accordance  with  the  direction  and  magnitude  of  each  force. 
In  this  case  the  effective  mass  acting  against  all  the  supports  should 
be  determined  by  adding  the  participation  factors  by  absolute  value 
(that  is.  ignoring  their  directions)  and  then  squaring  the  sum.  Notice 
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that  with  independent  motions  of  multiple  supports,  the  reactions 
occur  and  the  supports  are  moved  whenever  a  mode  is  accelerated, 
whether  the  acceleration  was  initiated  by  a  shock  delivered  through  a 
particular  support  or  not. 

Using  the  square  of  the  sum  of  the  absolute  values  of  the  mass- 
normalized  participation  factors  as  a  measure  of  the  mass  of  a  normal 
mode  of  vibration  is  the  only  substantial  change  that  the  RIDDAM 
makes  in  the  IDDAM.  The  IDDAM  distributed  the  mass  of  each  mode 
among  the  supports  in  accordance  with  the  shock  severity  at  each 
support.  In  the  RIDDAM  each  mode  has  a  mass  that  is  independent  of 
the  shocks  that  are  applied. 


In  the  original  DDAM.  the  shock  motion  of  a  ship  was  considered 
to  be  a  rigid-body  translation  in  a  single  direction.  All  the  supports 
for  an  item  of  equipment  underwent  the  same  motion  and  the  masses 
of  the  modes  were  the  effective  masses  reacting  against  that  common 
motion.  In  the  RIDDAM  each  support  is  taken  as  attached  to  a 
structure  that  can  move  independently  of  the  other  supports  and  the 
effective  mass  of  each  mode  reacts  against  all  of  the  supports  in 
whatever  direction  applies  at  each  support. 

It  is  recognized,  of  course,  that  points  located  near  one  another 
on  the  structure  of  a  ship  are  neither  lied  rigidly  together  l DDAM)  nor 
completely  independent  (RIDDAM).  The  object  of  the  procedures 
adopted  in  each  case  is  to  provide  a  reasonable  and  simple  method  of 
estimating  the  effects  of  the  reaction  forces  without  having  to  extend 
the  analysis  into  the  details  of  the  supporting  structure. 


The  example  shown  in  Equations  10  to  M  was  for  a  single  mass 
moving  in  a  single  direction.  Later  examples  will  show  that  the 
procedure  of  adding  influence  coefficients  by  absolute  value  and 
squaring  the  sum  gives  an  effective  mass  for  a  mode  even  when  a 
structural  model  includes  translations  in  different  directions  or 
rotations  about  different  axes. 
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The  IDDAM  inputs  represent  linear  accelerations  as  a  function 
of  effective  mass  and  frequency.  Applied  through  different  supports, 
these  accelerations  can  induce  rotational  and  cross-axis  responses  of  a 
structure.  However,  no  inputs  are  supplied  to  show  how  the  response 
of  a  structure  to  a  rotation  of  one  of  its  supports  might  depend  on  the 
mass  moment  of  inertia  about  that  support. 

The  RIDDAM  uses  the  inputs  developed  for  the  IDDAM.  As  a 
result  of  the  lack  of  rotational  inputs,  it  is  necessary  that  the  structure 
be  modeled  so  that  its  response  depends  on  translational  inputs  only 
This  can  be  done  in  either  of  two  ways:  A  support  can  be  taken  as 
irroutional,  so  that  moments  applied  to  it  do  not  cause  any  rotation 
Alternately,  the  joint  at  the  support  can  be  taken  as  pinned,  so  that 
the  structure  can  rotate  freely  about  that  point  without  producing  any 
rotation  of  the  supporting  structure 


INPUTS 


Differential  equations 


The  inputs  to  the  analysis  consist  of  predetermined  solutions  to 
the  equations 


Me  *  M  P2  z  »  M  P2  u, 
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where  dt)  is  the  acceleration  and  zlt)  is  the  displacement  of  an 
oscillator  having  mass  M  and  frequency  P  Here  uU'i  is  the 
displacement  of  the  support  for  the  oscillator,  including  both  the 
motion  induced  by  an  external  force  and  the  motion  resulting  from  the 
reaction  of  the  mass  of  the  oscillator  against  its  support. 

Because  of  the  reaction,  the  acceleration  of  the  oscillator  varies 
with  its  mass  as  well  as  with  its  frequence  In  the  DDAM.  all  the 
support  motions  are  taken  as  identical,  so  onlv  one  function  is  needed 
for  the  base  motion  utt),  but  in  the  IDDAM  the  supports  can  move 
independently  and  there  is  a  separate  Equation  1.5  for  the  effect  of 
each  support  on  each  mode 


The  equation  for  the  response  of  one  of  the  normal  modes. 


b  ♦  P2  y  -  P2  rT  u,  (16) 

uses  a  row  from  the  influence-coefficient  matrix  to  take  a  linear 
combination  over  the  motions  of  ail  of  the  supports  and  uses  that 
combination  to  find  the  response  of  the  mode.  The  process  can  be 
reversed  by  calculating  responses  first  and  then  taking  the  linear 
combination  to  obtain 

b  -  rT  c.  1 1 7 1 

where  b(t)  is  the  acceleration  of  the  mode  and  c(M,  P.  ti  is  a  column 
vector  of  the  solutions  to  Equation  15  for  the  mass  and  frequency  of 
the  mode  and  for  the  shock  applied  to  each  of  the  supports 

If  the  response  accelerations  c(M,  P,  t)  for  each  support  are 
already  provided.  Equation  17  will  give  the  responses  for  the  modes  of 
the  substructure  without  the  necessity  of  solving  any  differential 
equations 

ACCELERATIONS.  FORCES.  STRESSES.  AND  DEFORMATIONS 
Accelerations 

Scale  the  matrix  of  raodeshapes  to 

A  *  Q  B,  I  1 8  1 

where  B(t)  is  a  diagonal  matrix  of  the  accelerations  bit  I  of  the 
individual  modes  and  A*u  is  now.  by  Equation  5.  a  matrix  whose 
columns  represent  the  accelerations  of  the  coordinates  of  the  structure 
in  each  mode.  It  is  convenient  at  this  point  to  keep  the  responses  in 
the  different  modes  separate  because  of  approximations  that  will  be 
made  later  in  combining  responses  in  different  modes. 


Multiply  the  matrix  of  accelerations  by  the  mass  matrix  M  to 


obtain 


G  -  M  A  (19) 

as  a  matrix  whose  columns  represent  the  forces  applied  to  the 
structural  nodes  by  the  accelerations  of  each  mode 

Stresses  and  deformations 

The  forces  will  produce  time-varying  stresses  and  deformations 
in  the  structure  that  are  proportional  to  the  combination  of  forces  m 
each  mode.  The  stresses  and  deformations  at  critical  points  on  the 
structure  can  be  used,  to  estimate  the  ability  of  the  structure  to 
survive  the  accelerations  given  by  the  inputs 

WARPAGE 

Additional  deformations 

The  displacement-balance  equation  (Equation  2  >. 

i  -  K  1  F  u  t'l  M  a  120) 

shows  that  the  response  xm  of  the  structure  depends  on  the  motions 
of  us  supports  uti)  as  well  as  on  us  acceleration  am  The  preceding 
modal  analysis  accounts  for  the  forces,  stresses,  and  deformations 
produced  by  acceleration  only  It  is  necessary  to  supplement  the 
acceleration-generated  responses  with  a  check  of  anv  additional 
deformations  (and  resulting  stresses)  that  mav  be  produced  by 
differences  in  motions  of  individual  supports 

Note  that  if  the  motions  of  the  supports  are  all  identical  'as 
assumed  m  the  DDAM>  there  will  a  displacement  response  but  no 
additional  deformation  of  the  structure  There  also  is  no  warping 
deformation  if  the  supports  are  not  redundant  or  if  the  motions  of  the 


supports  can  be  described  by  a  combination  of  rigid -body 
displacements  and  rotations. 


The  structure  will  resist  warping  deformation  by  generating 
static  forces  that  tend  to  decrease  the  differences  between  the  motions 
of  its  supports.  Again,  as  in  the  case  of  the  effect  of  modal  mass  on 
acceleration  responses,  an  estimate  of  the  effect  of  the  reaction  forces 
can  be  made  in  terms  of  an  experimental  or  analytical  decrease  in 
warping  with  increasing  static  reactions  of  a  structure  against  relative 
motions  of  a  standard  or  typical  supporting  structure 

SUMMARY 

The  preceding  theory  shows  how  the  responses  of  substructures 
can  be  estimated  under  standard  conditions  while  mounted  to  a 
particular  supporting  structure.  The  characteristics  of  the  supporting 
structure  are  evaluated  by  determining  the  responses  it  can  produce 
in  oscillators  having  different  masses  and  frequencies  and  the  static 
warpage  it  can  produce  in  redundant  supports  with  different 
reactions 

The  structure  of  interest  is  then  analyzed  into  oscillators  having 
different  masses  and  frequencies  and  into  reactions  against  relative 
motions  of  us  supports  Direct  comparison  with  the  responses 
determined  for  the  supporting  structure  then  allows  responses  of  the 
structure  to  be  estimated 

The  theory  has  the  potential  for  being  exact  if  the  supporting 
structure  were  identical  from  one  substructure  to  the  next.  In 
application,  however,  supporting  structures  tend  to  be  similar  but  not 
identical.  Then  the  accuracy  of  the  method  depends  on  how  much 
variety  is  allowed  in  supporting  structures  before  a  new  set  of 
characteristics  must  be  applied 
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APPLICATION 


UNITS 

Convenience 

The  IDDAM  is  based  on  user-friendly  engineering  units  in  which 
Newton  s  law  is  satisfied  for  forces  in  pounds,  masses  in  pounds,  and 
accelerations  in  multiples  of  the  acceleration  of  gravity.  The 
acceleration  of  gravity  is  defined  as  386  inches  per  second  per  second 
(980-4-4  meters  per  second  per  second  with  one  inch  equal  to  0.0254 
meters).  A  pound  of  force  is  then  4.447  Newtons,  while  a  pound  of 
mass  is  0  4536  kilograms. 

Inputs 

The  tabulated  inputs  for  the  IDDAM  give  accelerations  m 
multiples  of  the  acceleration  of  gravitv  ig)  as  a  function  of  weights  in 
pounds  and  frequencies  in  Hertz  i cycles  per  second). 


Most  computer  programs  will  not  allow  forces  and  masses  both 
to  be  entered  in  pounds  A  common  convention  for  using  such 
programs  m  the  inch-pound-second  svstem  is  to  enter  masses  in  units 
of  386  pounds  That  is.  weights  in  pounds  are  divided  by  386  to 
obtain  masses  for  entry  into  the  computer  program  Forces  can  then 
be  entered  in  pounds,  with  stiffnesses  in  pounds  per  inch  and  elastic 
moduli  in  pounds  per  square  inch.  There  is  no  accepted  name  for  the 
unit  of  mass  that  weighs  386  pounds,  although  the  term  siinch  has 
been  suggested  as  a  parallel  to  the  term  slug'  that  is  used  for  the 
similar  unit  m  the  foot-pound-second  svstem 


SCALING  OF  MODESHAPES 


Some  computer  programs  apply  scale  factors  to  the  modeshapes 
to  produce  a  modal  matrii  QS,  where  S  is  a  diagonal  matrix  of  scale 
factors.  One  popular  method  of  scaling  makes  the  largest  element  in 
each  column  of  QS  have  the  value  unity,  for  example.  The  theory 
shown  here  was  based  on  mass-normalized  modeshapes  in  which 
qtmq  -  I  was  a  unit  matrix.  Modifications  are  needed  if  modeshapes 
have  been  scaled  otherwise. 


With  scaled  modeshapes, 
(QS)T  M  QS  -  S2 


is  diagonal  but  not  a  unit  matrix.  Also, 


(QS|T  K  QS  «  S2  p2  122) 

is  diagonal  but  its  elements  are  no  longer  the  squares  of  the  natural 
frequencies.  Using  the  scaled  modeshapes  in  the  transformations 


a  -  QS  b 


x  -  QSy 


converts  the  force-balance  equation 


M  a  '  Hi  -  F  u 


S2  b  *  S2  P2y  -  (QS)T  F  u 
as  the  analog  of  the  normal-mode  equations. 


net 


The  elements  of  S2  in  Equation  25  have  been  referred  to  as 
masses  of  the  normal  modes  and  even  used  as  if  they  were  masses  in 
some  calculations.  One  popular  computer  program  refers  to  them  as 
"generalized  masses'"  in  its  printed  output.  The  elements  of  S2  have 
completely  arbitrary  positive  values  that  depend  only  on  what 
criterion  was  used  to  scale  the  modeshapes  and  have  no  physical 
significance 


The  scale  factors  contained  in  S  must  simply  be  removed  step- 
by-step  as  the  analysis  proceeds.  First,  divide  by  the  elements  of  S2 
to  obtain  the  squares  of  the  natural  frequencies  as 

p2  .  s-2(S2P2).  (26) 

Next,  recognize  two  kinds  of  participation  factors.  The  mass- 
normalized  participation  factor  must  be  recovered  from 

Rl  -  S-MP-2  (0S)tF],  (27) 

The  responses  of  the  scaled  modes,  however,  must  be  found  from  the 
equation 

b  ♦  p2y  ,  p2  [P'2  S-2  (QS)tF1  u,  1 28) 

so  that  the  participation  factor  for  this  operation  must  be  obtained 
from 

r,T  ,  p-2s-2  (QS)t  F.  <2 9) 

Notice  that  any  computer  program  that  displays  natural 
frequencies  must  have  calculated  a  value  of  S2  in  order  to  evaluate 
Equation  26,  so  it  should  not  be  necessary  to  recalculate  the  scale 
factors  (Equation  2 1 )  in  order  to  find  what  factors  to  remove. 
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Many  computer  programs  begin  a  normal-mode  analysis  by 
choosing  a  modeshape  matrix  that  reduces  qtmq  to  a  unit  matrix 
With  such  programs  the  mass-normalized  modeshapes  are 
fundamental  and  subsequent  scaling  may  be  an  option  that  need  not 
be  selected. 

If  only  scaled  modeshapes  are  available  from  a  particular 
program,  it  may  be  more  efficient  to  remove  the  scale  factors  directly 
from  the  modeshapes  by  dividing  each  column  of  QS  by  its  element  of 
S  and  then  proceeding  in  a  mass-normalized  fashion,  rather  than 
dividing  repeatedly  by  the  scale  factors  as  results  are  being  calculated 

SUPPRESSION  OF  TIME  HISTORIES 

Complication 

The  theory  indicates  that  a  supporting  structure  must  be 
characterized  by  a  set  of  time-history  accelerations  ciM.  P.  u  that 
differ  for  each  mass  M  and  frequency  P  of  oscillator  that  may  be 
supported  from  it.  Moreover,  the  time  histories  cannot  be  oscillations 
at  the  frequency  P,  because  P  is  a  fiied-base  frequence  and  the 
essential  feature  of  the  supporting  structure  is  that  it  cannot  be 
considered  to  be  a  fixed  base. 

In  addition,  the  response  of  a  multimode  structure  with 
multiple  independent  supports  must  be  determined  bv  superposing 
separate  time-historv  responses  for  each  mode  and  each  support 
Such  combinations  were  impractical  for  routine  analyses  when  DDAM 
was  under  development  about  1960. 

Simplification 


The  IDDAM  follows  the  lead  established  by  the  DDAM  in 
suppressing  the  time  variation  of  the  response  acceleration  it 
characterizes  supporting  structures  in  terms  of  the  peak  'largest 
absolute  value)  of  acceleration  that  they  can  produce  for  an  oscillator 
of  mass  M  and  fixed-base  frequency  P 


Suppressing  the  time  histories  allows  supporting  structures  to 
be  described  by  double-entry  tables  showing  peak  acceleration  as  a 
function  of  mass  and  frequency,  and  eliminates  the  superposition  of 
time  histories  for  multiple  modes  and  multiple  supports  in 
determining  responses.  It  leads  to  problems,  however,  in  estimating 
the  combined  effects  of  several  modes  whose  individual  responses  are 
given  only  in  terms  of  their  peak  values. 

TIME  HISTORIES  AND  COMBINATION  RULES 


Uooer  and  lower  bounds 

It  is  clear  that  if  the  peak  responses  of  several  modes  were  all 
to  occur  in  the  same  direction  at  the  same  instant  of  time,  the 
combined  response  would  be  the  sum  of  the  individual  peaks  and  that 
this  sum  would  be  an  upper  bound  for  a  time-history  combination  of 
responses.  A  lower  bound  may  also  exist  if  one  of  the  peak  values  is 
large  enough  to  dominate  the  combination  even  if  all  the  other  peaks 
happened  to  occur  simultaneously  in  the  opposite  direction.  Either 
occurrence  is  possible  but  unlikely  in  a  sum  of  time  histories. 

Combination  rules 

Simple  and  arbitrary  rules  can  be  defined  to  estimate  a  peak 
from  a  combination  of  responses  that  are  described  only  by  their 
individual  peak  values.  Any  rule  adopted  ought  to  define  a  combined 
peak  falling  in  the  range  from  the  lower  to  the  upper  bound.  The  rule 
should  select  a  likely  value  m  that  range  and  should  be  adjusted  to 
avoid  too  much  overconservatism  or  underconservatism. 

The  exact  form  of  a  combination  rule  is  not  important  in  areas 
of  a  structure  where  one  mode  predominates,  so  that  the  upper  and 
lower  bounds  do  not  differ  greatly  A  combination  of  engineering 
judgment,  common  sense,  and  experimental  data  are  needed  to 
determine  combination  rules  when  several  modes  have  comparable 
peak  values  so  that  there  is  a  large  spread  between  upper  and  lower 
bounds.  Under  this  condition  the  best  rule  would  be  one  that  avoids 
obvious  conditions  of  great  overestimates  or  underestimates,  rather 


than  a  rule  that  attempts  to  recreate  an  actual  combination  of  time- 
history  responses. 


COMBINING  PEAK  RESPONSES  FROM  DIFFERENT  SUPPORTS 
Rule 

A  suitable  rule  for  estimating  the  peak  acceleration  b  of  one 
mode  from  a  collection  of  peak  response  accelerations  c  at  different 
supports  is 

b  *  rT  c.  1 30  i 

where  is  the  row  for  the  mode  from  the  participation-factor  matrix 
Discussion 

Equation  30  is  identical  to  Equation  17  that  was  used  for 
combining  time  histories.  It  presumes  that  the  peak  responses  of  the 
mode  produced  by  the  motions  of  each  support  all  occur  at  the  same 
instant  of  time  and  in  directions  given  by  the  signs  of  the  elements  of 
the  participation  factor  Such  a  combination  could  actually  occur 
under  torpedo  impact  if  the  torpedo  were  to  strike  at  the  center  of  a 
symmetric  array  of  supports,  producing  large  responses  to  symmetric 
modes  of  the  structure  and  no  response  for  antisymmetric  modes 

Equation  30  also  reconciles  the  RIDDAM  with  the  DDA.V1  by 
matching  its  results  with  those  of  the  DDAM  for  conditions  in  which  all 
of  the  supports  would  have  identical  motions  in  a  single  direction. 

In  other  cases  the  rule  provides  an  estimated  peak  that  falls 
between  the  upper  and  lower  bounds  for  the  combination. 

COMBINING  PEAK  RESPONSES  FOR  MODES  OF  DIFFERENT 
FREQUENCIES 

The  NRL  Sum 

An  estimate  of  the  peak  stress  or  deformation  from  a 
combination  of  peak  stresses  or  deformations  from  modes  of  different 
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frequencies  is  obtained  by  electing  the  response  having  the  largest 
peak  in  absolute  value  and  adding  to  it  the  square  root  of  the  sum  of 
the  squares  of  the  peak  values  of  the  other  responses.  This  method  of 
combination  is  usually  called  the  NRL  Sum. 

Estimated  peaks  obtained  from  the  NRL  Sum  are  always 
positive  numbers,  but  are  to  be  interpreted  as  representing  peaks  that 
may  occur  in  either  direction. 

Warning 

The  NRL  Sum  is  a  nonlinear  combination  and  must  always  be 
the  last  step  in  any  calculation.  In  particular  if  the  NRL  Sum  is  used 
to  estimate  a  peak  acceleration  or  a  peak  force  from  a  combination  of 
modal  peaks,  that  estimated  peak  acceleration  or  force  cannot  be  used 
subsequently  to  find  stresses  or  deformations 


stimati 


)eak  stresses  and  deformations 


Equations  18  and  19  were  especially  written  in  a  form  to  keep 
the  responses  of  the  structure  in  each  of  its  modes  separate  down  to 
the  point  of  calculating  stresses  and  deformations.  If  B  is  a  diagonal 
matrix  of  the  peak  accelerations  b  in  each  mode,  the  peak 
accelerations  of  the  structure  in  each  mode  are  obtained  by  scaling  the 
modeshape  matrix  according  to 


The  peak  forces  on  the  structure  in  each  mode  are  given  by  the 
columns  of 

G  -  M  A  (32' 

The  stresses  and  deformations  for  each  mode  are  obtained  by  solving 
a  separate  static  prohlem  for  the  peak  loads  g  that  appear  in  each 
column  of  the  matrix  of  forces  The  loads  mav  be  positive  or  negative, 
depending  on  the  signs  of  the  participation  factors  and  the 
modeshapes;  these  signs  must  be  preserved  through  Equations  31  and 
32 
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The  strain  energy  associated  with  a  peak  response  is 
proportional  to  the  square  of  the  strain  multiplied  by  the  elastic 
modulus  of  the  material,  or  proportional  to  the  square  of  the  stress 
divided  by  the  elastic  modulus.  Summing  the  squares  of  the  peak 
responses  thus  gives  a  number  proportional  to  the  total  energy  a 
particular  part  of  the  structure  must  accept  The  square  root  of  that 
sum  represents  the  single  deformation  or  stress  that  would  account 
for  the  same  energy 

The  NRL  Sum  provides  an  estimated  peak  response  that  is 
larger  than  the  energy -equivalent  peak,  but  not  larger  than  the  upper 
bound 

SELECTION  RULES  FOR  THE  NRL  SUM 


Peaks  from  the  following  sources  are  combined  into  a  single 
NRL  Sum: 

1  Stresses  and  deformations  from  normal  operation  of 
equipment 

2.  The  warping  deformations  and  stresses  from  differences  in 
motions  of  multiple  supports 

3  Stresses  and  deformations  from  some  but  not  usually  all  of 
the  normal  modes.  The  NRL  sura  overestimates  the  energies  lor  the 
modes  that  are  included:  this  overestimate  can  be  considered  as  an 
allowance  for  additional  modes  of  the  structure  that  were  not  included 
in  the  sum 


1  For  simple  structures  that  can  be  modeled  with  six  or  fewer 
degrees  of  freedom,  the  three  modes  of  lowest  frequency  are  usually 
sufficient  for  determining  an  NRL  Sum. 

2.  Structures  of  moderate  complexity  'up  to  bO  degrees  of 
freedom)  need  only  the  lower -frequency  modes  to  be  included  in  the 
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NRL  Sum.  About  half  as  many  modes  as  degrees  of  freedom  is  an 
appropriate  combination  in  most  cases. 

3.  Complicated  models  with  very  many  degrees  of  freedom  (as 
may  be  obtained  from  a  finite-element  model)  may  generate  a  “modal 
thicket"  in  which  many  of  the  low-frequency  modes  represent 
vibrations  of  lightweight  parts  of  the  structure.  Here  it  may  be  better 
to  select  modes  in  accordance  with  modal  weight  rather  than 
frequency.  Modes  having  weights  at  least  2  percent  of  the  weight  of 
the  complete  model  are  most  likely  to  be  significant. 

Prestresses  and  bolted  joints 

Built-in  stresses,  including  prestresses  in  bolts,  are  not  included 
in  the  NRL  Sum  of  stresses.  Bolted  joints  are  expected  to  respond  to 
tensile  'loads  by  decreasing  the  clamping  force  rather  than  by 
stretching  the  bolts.  Special  checks  need  to  be  made,  however,  to  be 
sure  than  the  clamping  force  is  not  exceeded  by  a  peak  load  or  that 
shear  stress  from  a  transverse  load  does  not  combine  with  the  tensile 
stress  in  a  bolt  to  produce  an  excessive  value  of  maximum  normal 
stress. 

SUMMARY 

The  procedures  for  applying  the  theory  to  practical  calculations 
with  the  RIDDAM  can  be  summarized  in  eleven  steps 

1.  Make  a  mathematical  model  of  the  structure. 

2.  Find  the  fixed-base  normal  modes  and  the  mass-normalized 
participation  factors  for  the  supports  of  the  model. 

3.  Add  the  participation  factors  for  each  mode  by  absolute 
value  and  square  the  sum  to  find  the  effective  mass  of  each  mode 

4.  Refer  to  tables  and  formulas  showing  peak  accelerations  as  a 
function  of  the  weight,  frequency,  and  the  location  of  an  oscillator 
relative  to  the  point  of  impact  of  a  torpedo 

V  Combine  the  product  of  the  participation  factors  and  the 
tabulated  accelerations  for  each  support  to  find  the  peak  acceleration 
of  each  mode 

6.  Scale  the  modeshapes  by  the  peak  accelerations  to  find  the 
peak  acceleration  of  each  point  on  the  structure  in  each  mode 


7.  Multiply  the  peak  accelerations  of  the  structure  by  the  mass 
matrix  to  find  peak  forces  in  each  mode. 

8.  Solve  the  static  problem  for  each  mode  to  find  peak  stresses 
and  deformations  at  critical  parts  of  the  structure. 

9.  Calculate  stresses  and  deformations  from  warping  of  the 
structure  caused  by  different  motions  of  redundant  supports.. 

10.  Use  the  NRL  Sum  to  combine  stresses  and  deformations 
from  operation,  warping,  and  selected  modes  at  critical  points  in  the 
structure. 

11.  Compare  the  estimated  peak  stresses  and  deformations 
with  allowable  values. 


SPECIAL  PROBLEMS 


DESIGN  SPECTRA 

The  accelerations  c(M.  P)  representing  the  peak  accelerations 
for  oscillators  having  masses  M  and  fixed-base  frequencies  P  can  be 
called  a  design  spectrum.  Each  point  in  a  design  spectrum  is  taken 
from  a  particular  time-history  response  of  a  massive  oscillator.  That 
time  history  has  been  suppressed  for  simplicity. 

The  special  problems  described  here  arise  from  the  suppression 
of  the  time  histories  associated  with  the  design  spectra. 

REPEATED  FREQUENCIES 

Time -history  combinations 
Suppose 

ci  *  Pj2  z\  -  P}2  u  (53 1 

and 

C2  *  P22  Z2  *  P22  u  (34) 

are  the  responses  of  two  oscillators  to  a  support  motion  u«u  If  the 
two  frequencies  are  the  same,  Pj  «  P?  -  P  the  time  histories  will  be 
identical  and  the  peak  accelerations  will  occur  simultaneously 

As  the  accelerations  are  scaled  by  participation  factor  and  by 
modeshape  and  then  used  to  calculate  forces,  stresses,  and 
deformations,  the  peaks  will  remain  coincident  in  time  but  may 
represent  responses  of  different  magnitude  occurring  in  either  the 
same  or  opposite  directions  The  combined  peak  is  given  by  the 
algebraic  sum  of  the  individual  peak  values 
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If  Modes  1  and  2  are  the  only  two  inodes  of  the  structure,  the 
NRL  Sum  will  add  their  peak  stresses  and  deformations  by  absolute 
value  (absolute  value  of  the  larger  plus  square  root  of  the  square  of 
the  smaller).  This  is  an  appropriate  combination  if  both  responses  are 
in  the  same  direction,  but  it  may  be  a  very  large  overestimate  if  the 
modes  actually  produce  responses  in  the  opposite  directions. 

If  the  structure  has  more  than  two  modes,  either  or  both  of 
Modes  1  and  2  may  be  part  of  the  square  root  of  the  sum  of  the 
squares  in  the  NRL  Sum.  Here  their  combined  response  will  be 
underestimated  if  thev  respond  in  the  same  direction,  and  will  be 
overestimated  if  they  respond  in  opposite  directions. 


Modes  with  identical  frequencies  should  have  their  peak 
responses  added  algebraically  rather  than  by  NRL  Sum 


CLOSELY-SPACED  MODES 


Formal  solutions  for  Equations  33  and  34  can  be  written 


Cl  =  Pi2  |  vl. t  )  COSlPl't  -  t  I)  dl 
0 


C2  =  P2^  J  vi t  i  cos(P2«t  -  l  »j  dt 
0 


where 


is  the  velocity  of  the  support. 


Time-history  combinations 

When  the  two  accelerations  are  scaled  by  participation  factor 
modeshape,  mass,  and  stress  factor  they  will  produce  a  contribution 

t 

sit)  -  Dj  \  v(t  *  cosfPiit  -  t  >j  dt 
0 


t 

*  D2  |  vi t  )  cos(P2<t  -  t  )]  dt  '38  1 

0 

to  a  particular  stress,  where  Dj  and  D2  are  combinations  of  all  the 
relevant  factors  for  the  two  modes. 

If  the  frequencies  are  Pi  -  P  -  E  and  P2  «  P  *  E.  with  an  average 
frequency  P  and  a  difference  of  2E.  the  stress  can  be  written 

t 

s'tl  *  1 D 1  *  D2 1  |  v(t  )  cos(E( t -t  1]  coslP't-t 'I  dt 
0 

t 

*  'Dj  -  D21  J  v(t  )  smiEu-t  1  sinfPlt-t  i|  dt  '  ><■»( 

0 


Approximation  for  closely-spaced  frequencies 

If  Et  is  less  than  tt/2.  the  factors  cos[E't-t  )J  and  sinlE't  1  i|  can 
be  treated  as  positive-valued  windows  applied  to  the  velocity  vt  1 


with  average  values  — ^ —  and  - ^ -  over  the  range  of 

integration' from  0  to  t.  The  averages  can  be  considered  as  weighting 
factors  applied  to  the  velocity  of  the  support.  The  first  average  begins 
with  a  value  of  unity  at  the  beginning  of  the  response  and  decreases 
with  increasing  time,  while  the  second  begins  at  zero  and  increases. 

A  criterion 

An  Et  of  tt/A  produces  an  average  of  0  900  for  the  first  window 
and  an  average  of  0.373  for  the  second.  For  frequency  differences  and 
times  up  to  this  limit  the  stress  mav  be  approximated  bv 

t 

sill  -  (Di  ♦  D2(  j  vi t  *  cos{P< t-t  f|  dt .  ( -40 1 

0 

to  reasonable  accuracv  That  is.  both  time  histories  approximate  the 
response  of  an  oscillator  at  the  average  frequency  P  and  their 
combined  stress  is  obtained  from  an  algebraic  sum  of  the  factors  D[ 
and  E>2.  as  was  the  case  for  identical  frequencies 

If  the  frequency  difference  2 E  is  3  Hertz,  for  example.  Equation 
40  can  be  used  to  estimate  combined  peak  stresses  for  times  up  to  SO 
milliseconds  alter  a  transient  motion  of  the  support  begins 

Further  approximations 

Convenience  is  served  il  the  responses  of  modes  having 
identical  or  closelv-spaced  frequencies  are  calculated  independent^'  m 
terms  of  their  weights  and  frequencies  Then  when  peak  responses 
are  to  be  combined  for  stresses  or  deformations  peaks  lor  pairs  or 
clusters  of  modes  having  near:\  the  same  frequencv  should  be 
combined  bv  algebraic  sure  u.  produce  a  single  response  That 
response  should  be  treated  as  the  response  of  a  single  mode  and 
combined  with  other  responses  bv  'he  NRL  Sum 

Measurements  and  analvses  >>i  the  responses  ol  oscillators  to 
torpedo  impact  against  a  submarine  regularlv  show  that  peak 


responses  occur  less  than  Si)  milliseconds  alder  the  initial  impact  The 


peak  stresses  and  deformations  from  pairs  or  clusters  of  modes  that 
fall  within  a  5-Hertz  range  can  then  be  combined  algebraically  rather 
than  by  NRL  Sum. 

RULES  OF  THUMB 

Suppression  of  the  time  histories  is  seen  to  require  a  number  of 
fairly-arbitrary  rules  concerning  selection  of  modes  and  methods  of 
estimating  peak  values  from  combinations  of  responses 

The  combination  rules  suggested  here  appear  to  be  based  on 
rational  criteria.  The  rules  were  chosen  mamiv  to  avoid  obviously- 
mappropriate  results  and  are  subject  to  revision  and  reinterpretation 
in  light  of  statistical  results  from  experimental  data. 
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EXAMPLES 


A  TEXTBOOK.  EXAMPLE 

Structure 

Rigid  masses  weighing  325  pounds  and  200  pounds  are 
supported  by  springs  having  stiffnesses  of  10,000  and  30,000  pounds 
per  inch,  as  shown  in  Figure  i.  Each  mass  is  prevented  from  rotating 
by  a  set  of  roller  guides  that  are  not  considered  part  of  the  supporting 
structure.  The  masses  are  connected  by  a  thm  and  massless  beam  30 
inches  long.  The  cross-section  of  the  beam  has  an  area  moment  of 
inertia  of  1.25  in4  and  the  beam  is  made  of  a  material  with  elastic 
modulus  30  million  pounds  per  square  inch. 

The  structure  is  devised  to  illustrate  the  calculations;  it  is  not 
intended  to  represent  anything  realistic. 

M35S  matrix 

For  coordinates  x j(t)  and  xj'ti  representing  displacements  of 
the  masses,  the  mass  matrix  is 


M 


T325  0-, 

L  0  200 J  P°unds- 


Ml' 


Stiffness  matrix 

Moving  xi  by  one  inch  and  holding  all  the  other  coordinates 
fixed  requires  a  force  of  10,000  pounds  to  stretch  the  spring  and  an 
additional  force  of 


12  (30E61  1  25 
(30)3 


1 6,666.67  pounds 


I42< 


to  bend  the  beam  as  a  double  cantilever  An  equal  and  opposite  force 
must  be  applied  to  X2  to  keep  the  other  end  of  the  beam  from  moving 


1 fc*l 


10,000 


325 


200 


30,000 


:igure  1.  Textbook  structure.  Rigid  masses  are  labeled 
in  pounas.  Massless  springs  are  labeled  in  pounds  per 
inch.  The  masses  are  connected  by  a  massless  thin  beam 
30  inches  long  with  area  moment  l  25  inches  to  the 
fourth  power  and  elastic  modulus  30  million  pounds  per 
square  inch.  The  masses  are  on  frictionless  guides  that 
prevent  them  from  rotating 
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Here  and  in  the  following  the  usual  computer  notation  of  E6  is  used, 
for  example,  to  indicate  multiplication  by  10  to  the  6th  power 

Moving  12  while  the  other  coordinates  are  held  fixed  stretches 
its  spring  and  also  bends  the  beam.  The  forces  from  unit  displacement 
of  one  coordinate  at  a  time  can  be  assembled  into  the  stiffness  matrix 


f  26666.67 
L  - 1 6666  67 


- 1 6666.671 

46666.67  J  Pounds  per  inch. 


A  one-inch  displacement  of  ui  (with  everything  else  held  fixed) 
will  put  a  force  of  10,000  pounds  onto  xi.  A  similar  displacement  of 
U2  will  put  a  force  of  30.000  pounds  onto  X2.  The  support-stiffness 
matrix  becomes 


0000 

0 


O' 

30000J 


pounds  per  inch. 


The  three  matrices  M.  K.  and  F  allow  the  structure  to  be 
described  by  the  force-balance  equation  Ma  •  Kx  =  Fu.  shown  earlier 
as  Equation  1  in  the  section  on  Theory. 


The  next  step  in  the  analysis  is  to  find  a  matrix  0  that  will  make 
QtMQ  be  a  unit  matrix  and  also  make  QTKQ  be  a  diagonal  matrix,  as 
shown  in  Equations  3  and  4  This  can  be  done  most  eas.ly  by  using  an 
available  computer  program,  such  as  in  the  example  shown  as  Case  1 
in  Appendix  B,  For  the  simple  structure  here,  the  modeshape  matrix  Q 
can  be  calculated  bv  hand,  using  the  formula  listed  in  Line  4  1  20  of  the 
computer  program  to  normalize  on  the  masses  and  the  formulas 
beginning  on  Line  3610  to  diagonalize  the  stiffnesses  The  modeshape 
matrix 


r0  05  198516  -0.0 1935 1 1 4"; 

^  "  lO 02466795  0 06626833- 

is  from  a  10-digit  pocket  calculator.  A  direct  check  shows  that  it 

reduces 


QT  M  Q 


r  1  000  0.0001 
Lo.ooo  j.oooj 


to  the  unit  matrix  and 


QTKQ 


'57717  0.000" 

.  0.000  257.668J 


to  a  diagonal  matrix. 

The  matrix  0  can  be  used  to  represent  the  accelerations  and 
displacements  of  the  structure  in  terms  of  normal-mode  coordinates 
by  a  -  Qb  and  i  -  Qy  as  shown  in  Equation  5 

Participation  factor 


The  participation-factor  matrix  is 


P-2  QT  F 


r  9.007  1 2.82  1 1 

l  -0  75  1  7.7  16j 


-  rt 


The  participation  factor  lets  the  equation  for  the  normal-mode 
coordinates  be  written  in  the  standard  form  b  -  P2y  =  P2RTu 
(Equation  7).  Each  element  of  RTu  is  the  effective  displacement  of  the 
base  of  one  of  the  normal-mode  oscillators 


Add  along  the  rows  of  the  participation-factor  matrix  hv 
absolute  value  and  square  the  result  to  obtain 


(9  007  *  1 2.82  1  ^2  =  475  pounds 
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for  Mode  1  and 

(0.751  «■  7.7 1 6 12  -  72  pounds  150) 

for  Mode  2.  This  is  the  new  procedure  of  the  RIDDAM,  as  shown  by 
example  in  Equation  14. 


When  masses  are  in  pounds  and  stiffnesses  in  pounds  per  inch, 
frequencies  are  given  by  the  formulas 

~  V386  (57.717)  =  25.8  Hertz  i  5  1' 


86  (257.668)  -  50.2  Hertz 


for  Modes  1  and  2.  The  normal-mode  equations  for  motioniess 
supports,  b  -  P2y  >  0,  identify  ihe  diagonal  elements  of  P2  as  squares 
of  the  fixed-base  natural  frequencies  of  the  modes.  The  factor  of  586 
is  the  acceleration  of  gravity  in  inches  per  second  per  second,  as 
needed  in  the  inch-pound-second  system  when  pounds  are  used  for 
both  masses  and  forces 


Inputs 


The  weights  and  frequencies  of  the  modes  serve  as  entries  to 
tables  showing  the  peak  accelerations  of  the  modes.  As  an  example, 
suppose  that  the  tables  specify  a  peak  acceleration  of  26  times  the 
acceleration  of  gravity  ig)  for  a  weight  of  476  pounds  and  a  frequence 
of  23  8  Hertz,  and  an  acceleration  of  1  19  g  for  72  pounds  at  50  2  Hertz 
These  are  values  at  the  point  of  impact  for  a  torpedo  At  a  point  on 
the  hull  30  inches  away  from  the  point  of  impact  the  response 
accelerations  would  be  1  1  70  g  and  53  55  g.  These  values  represent 
peak  accelerations  of  oscillators  of  mass  M  and  fixed-base  frequenev  p 
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measured  during  tests  with  torpedo  impacts  or  calculated  from 
Equation  1 5  during  a  computer  simulation  of  an  impact. 

Peak  accelerations  of  the  modes 

If  the  torpedo  impacts  in  way  of  uj,  with  U2  located  30  inches 
away  from  the  impact,  the  peak  accelerations  of  the  modes  are 
obtained  by  multiplying  the  tabulated  responses  by  the  participation 
factors,  or 

9.007  (26)  ♦  12.821  (1 1.70)  =  38^.1 9  g  (33) 

for  Mode  1  and 

-0.751  (1  19)  *  7.716  (53.55)  =  323.82  g  (54) 

for  Mode  2. 

Equations  53  and  54  correspond  to  applying  the  participation 
factors  to  the  predetermined  responses,  rather  than  applvmg  the 
participation  factors  to  the  support  motions  and  then  calculating  the 
responses.  The  two  procedures  give  equivalent  results  for  time- 
history  calculations,  as  explained  in  Equations  16  and  17.  The  same 
combination  method  is  recommended  for  peak  values  in  Equation  30  if 
the  time  histories  have  been  suppressed. 

Peak  accelerations  of  the  trasses 

Scale  the  shape  for  the  first  mode  (the  first  column  of  Q)  by 
384.19  and  the  shape  for  the  second  mode  by  323.82  to  obtain 

I"  1 9.97  -6  27; 

A  =  L  9.48  21.45J  g 

as  the  peak  accelerations  of  the  masses  in  each  mode.  The  first 
column  of  A  shows  the  peak  accelerations  of  xj  and  x?  in  Mode  1  and 
the  second  column  shows  their  peak  responses  in  Mode  2.  Peak 
accelerations  are  showrn  separately  for  each  mode,  as  recommended  in 
Equation  3 1 . 


Multiply  the  masses  (pounds)  and  the  accelerations  tg)  to  find 


n 

ft 


„  .  I" 6490  -2038] 

A  -  i  4290  1  Pounc*s 


as  the  peak  forces  that  were  applied  to  each  mass  in  each  mode. 


Each  mode  is  treated  as  a  static  problem  to  check  for  stresses 
and  deformations  at  critical  points  on  the  structure.  Suppose  that  the 
critical  points  are  the  bending  moments  at  each  end  of  the  beam  and 
the  tensile  loads  applied  to  bolts  restraining  the  springs  A 
redundant-structure  analysis  provides  deflections 


r0.3-460 
[0  1642 


-0.0243  ' 


0.0832 


i  inches 


for  each  mass  in  each  mode  when  the  forces  of  Equation  56  are 
applied  as  static  forces  to  the  structure  The  deflections  produce 
forces  of  3460  pounds  at  ui  in  Mode  1  and  -243  pounds  in  Mode  2.  At 
U7  the  forces  are  4926  pounds  in  Mode  1  and  2496  pounds  in  Mode  2 


Relative  motions  between  the  supports  will  distort  the  structure 
through  its  redundant  supports.  Suppose,  for  example,  that  the 
warping  formulas  specify  a  displacement  of  0  200  inch  at  the  impact 
point  and  a  displacement  of  0  090  inch  a  distance  of  30  inches  from 
the  impact  point  A  redundant-structure  analysis  shows  that  x\  will 
have  an  equilibrium  displacement  of  0  1431  inch  and  x?  a 
displacement  of  0  1090  inch  under  these  support  displacements  The 
resulting  tensions  in  the  hold-down  bolts  are  -569  pounds  a:  Support 
1  and  an  equal  and  opposite  569  pounds  at  Support  2 


xv: ; 


Kv 
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Loads  on  the  supports  are  to  be  combined  by  NRL  Sum,  subject 
to  the  selection  rules  described  earlier  under  Application."  The  forces 
from  warping  and  two  modes  (assuming  no  operating  stresses) 
combine  by  the  rule  of  the  largest  plus  square  root  of  the  sum  of  the 
squares  in  the  form 


3460  -  V(- 243)2  -  (-569)2  .  4079  pounds  (38) 

as  an  estimate  of  the  peak  force  applied  to  the  bolts  at  uj.  and 

4926  *  \l( 2496)2  ♦  (569)2  ,  7486  pounds  (59) 

for  the  bolts  at  U2. 

Continuing 

The  peak  forces  of  Equations  58  and  59  are  interpreted  as 
applying  in  either  direction  to  the  hold-down  bolts.  The  bolts  need  to 
be  tight  enough  so  that  the  peak  force,  applied  in  tension,  will  not 
cause  the  joint  to  separate.  Similar  calculations  allow  peak  moments 
at  the  ends  of  the  beam  to  be  estimated  for  comparison  with  allowable 
moments. 

The  calculations  need  to  be  repeated  for  a  case  in  which  a 
torpedo  impacts  in  way  of  U2  rather  than  uj  The  most  severe 
conditions  can  be  taken  as  those  in  which  one  of  the  supports  for  an 
item  is  directly  inboard  of  the  point  of  impact  Then  a  check  of 
impacts  at  each  support  is  sufficient  to  evaluate  the  ability  of  the 
equipment  to  withstand  an  impact  at  any  point  near  it. 

AN  EXAMPLE  WITH  CROSS-AXIS  RESPONSE 

Structure 


A  point  mass  is  supported  by  two  equal  springs  battered  at 
angles  H  above  and  below  the  horizontal,  as  shown  in  Figure  2.  The 
mass  can  move  both  horizontally  ( x  1 )  and  vertically  ix?1  in  response  to 


the  horizontal  motions  of  its  two  supports.  This  is  another  textbook 
example  not  intended  to  represent  a  realistic  structure. 

In  this  and  in  the  following  two  examples  the  analysis  is 
presented  in  a  form  parallel  to  that  used  for  the  first  example,  but 
without  back-references  to  the  sections  of  the  report  on  theory  or 
application.  Please  refer  to  corresponding  sections  of  the  first  example 
for  appropriate  references. 

Matrices 

The  mass  matrix. 


M 


f  m  O' 

L  0  mJ  ' 


i 

i  6  0 )  « 


represents  motion  of  the  mass  in 
directions.  The  stiffness  matrix 

T  2  k  cos2/ H  i  0 

K  "  L '  0  2  k  sin2(H)  J  101 1 

includes  factors  for  the  components  of  force  and  deflection  in  each 
spring  produced  bv  small  displacements  of  the  mass  in  each  direction 
The  support-stiffness  matrix 


either  of  the  two  orthogonal 


l  L  l 


f  k  cos2< Hi 
L  -k  simHi  cosih'i 


k  cos2(  H  i 
k  sinlHi  cos/Hl 
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shows  the  horizontal  and  vertical  force  per  displacement  on  the  mass 
produced  by  small  horizontal  displacements  of  each  support. 

Normal  modes 


The  modeshape  matrix 


normalizes  the  mass  matrix  to  a  unit  matrix,  gives  a  frequencv- 
squared  matrix 


qtkq 


i *2  (k/m)  cos2(H)  0 


2  i k/m)  sin2(H)J 


and  a  participation-factor  matrix 


p-2  qT  p 


'/m/2 


V  m/2 


L  -(v m/2 ')  cotiH)  i'/m/2)  cotiH  l 


Mode  1  (horizontal  motion)  has  effective  mass  m,  as  obtained 
from  the  squared  sum  of  the  first  row  of  the  participation-factor 
matrix.  Mode  2  (vertical  motion'  has  effective  mass  m  cot^H1.  with 
the  square  of  the  cotangent  representing  the  effect  of  the  lever  arms 
that  act  to  magnify  or  diminish  the  effects  of  a  vertical  force  on  the 
horizontally -directed  supports  and  also  act  to  magnifv  or  dimmish  the 
motions  of  the  mass. 


Relative  horizontal  motions  between  the  two  supports  will 
simply  displace  the  mass  in  the  vertical  direction  without  producing 
any  stresses  or  deformations  (The  joints  at  the  supports  are  pinnec  ; 


The  simple  example  here  shows  that  the  matrix  formalism  is 
not  limned  to  unidirectional  responses  If  cross-axis  responses  are 
expected  from  unequal  motions  of  supports,  thev  should  be  included 
in  the  model  of  the  structure,  and  can  be  treated  on  the  same  basis  as 
on-axis  responses 


AN  EXAMPLE  WITH  ROTATIONAL  INERTIA 


Structure 

A  steel  plate  is  bridged  between  the  flanges  of  two  frames  to 
carry  a  weight  made  from  a  stack  of  steel  plates,  as  shown  in  Figure  5 
Four  structures  like  this  were  built,  installed  in  hulls,  and  tested  with 
torpedo  impacts.  The  one  analyzed  here  was  tested  in  1985  on  the 
SITV  (Submarine  Impact  Test  Vehicle)  There  it  served  as  simulated 
equipment  for  research  purposes 


The  structure  was  especially  designed  to  be  symmetric  and 
nonredundant  with  two  simple  modes  of  vibration  In  one  mode  the 
weight  would  translate  as  a  rigid  body  and  bend  the  plate  as  a  beam 
The  stiffness  of  the  plate  1 33  inches  iong.  22  inches  wide,  and  1  inch 
thick,  with  equal  loads  6.5  inches  from  each  end)  can  be  obtained  from 
handbook  formulas  as 


12  (28E6)  22  (1)3/ 12 
(6.5)2  (3(35)  -  4(6.5)! 


1 99 .7 24  pounds  per  inch 


The  stack  of  plates  15  59  pounds'  and  the  center  20  inches  ol  the 
flexible  plate  (124  pounds)  both  participate  in  the  bouncing  motion, 
while  the  6.5  inches  at  each  end  of  the  flexible  plate  connect  from  the 
fixed  support  to  the  moving  load  sav  one-third  of  81  pounds'  6 
rough  estimate  of  the  total  mass  m  the  bouncing  mode  is  thus  Mild 
pounds  The  frequency  should  be 


386  (199724' 
f  101 0 


44  nertz 


Analysis  of  measurements  from  the  SITV  test  showed  that  the 
bouncing  mode  had  a  fixed-base  frequency  of  44  Hertz 
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Figure  3  Structure  with  rotational  inert-a  Tre  weight 
(a  stack  of  steel  plates)  is  Doited  to  a  flexible  plate  that 
Dridges  Detween  the  flanges  of  two  successive  ring 
stiffeners  of  a  submarine  hull.  Half-round  spacers  (not 
illustrated)  provided  nearly  moment-free  connections  at 
the  ends  of  the  flexible  plate  and  at  the  connections 
between  the  plate  and  the  weight 


!*] 


If  the  weight  rotates  about  an  axis  through  the  center  of  the 
plate  and  parallel  to  its  ends,  each  half  of  the  plate  will  bend  as  a 
simply-supported  beam  with  a  point  load.  Handbook  formulas  give 
the  stiffness  of  each  half  as 

3  (28E6)  (22  (1)3/12]  16.5  . 

I ^  ^7  ( ]  ^  c  7  £  ^2 - *  601,420  pounds  per  inch.  (68; 


Two  such  loads  are  each  10  inches  from  the  center  of  the  plate,  and  a 
small  rotation  produces  deflections  of  10  inches  per  radian  The 
angular  stiffness  is  thus  120.284  million  pound-inches  per  radian 

The  mass  moment  of  inertia  for  rotations  about  an  axis  through 
the  center  of  the  plate  can  be  obtained  by  starting  with  the  moment 
about  the  center  of  gravity  of  the  weight  (a  block  25. 5  by  7  by  17 
inches), 


(25.5)2  *  (7)2 


-  50,055  pound-in2 


Transfer  to  a  parallel  axis  through  the  middle  plane  of  the  plate  i9 
inches  from  the  center  of  gravity  of  the  weight)  by 

850  (9)2  =  2  1 .475  pound-m2,  (70) 

and  add  an  estimate  for  the  center  part  of  the  plate 

.  „  (20  (2  -  ( 1 12 

124  7T  =  4,142  pound-m2.  (7p 


to  obtain  a  total  moment  of  inertia  of  75.682  pound-m2,  Frequencv  is 
then 


1  ,  ;'386  ( 120.284E61 


,25  Hertz. 


Data  from  the  test  showed  a  fixed-base  frequency  of  13?  Hertz  for  the 
rocking  mode 


Matrices 


The  mass  matrix, 


rioio 


M  '  l  0  75682J 

contains  a  mass  in  pounds  and  a  mass  moment  of  inertia  in  pound- 
inches  squared  The  stiffness  matrix. 


199724 


1 20.284E6.' 


has  elements  in  pounds  per  inch  and  in  pound -inches  per  radian.  The 
support-stiffness  matrix 

r  99862  9986 2; 

b  ”  L3.645E6  -3.645E6J 

shows  that  a  one-inch  displacement  of  both  supports  would  apply  a 
force  of  199,724  pounds  to  the  translational  coordinate. 
Displacements  of  - 1  and  -1  inch  would  cause  a  rotation  of  1/16  5 
radians  and  apply  a  moment  of  120  284E6/16  5  pound-inches  to  the 
rotational  coordinate 


The  modeshape  matrix 

[-3M6584E-3  0  ' 

Q  "  i  0  3  63499E-3-  76 

normalizes  the  mass  matrix  to  a  unit  matrix  and  gives  a  frequency- 
squared  matrix 


177) 


QT  K  Q  - 


f  197.746  0  000  1 

L  0.000  1589.334J 


=  P2 


that  gives  the  frequencies  of  44  and  125  Hertz. 

Participation  factor 

The  participation-factor  matrix  is 

d  2  at  c  "  1 5  8903  15.8903:  dT 

P  Q  F  “  i.  8.3365  -8.33652  ?  RT  ' 7 '  1 

Add  along  the  rows  by  absolute  value  and  square  the  sums  to  find 
that  Mode  1  (bouncing)  has  effective  weight  1010  pounds  and  Mode  2 
(rocking)  has  effective  weight  278  pounds 


From  this  point  onward,  the  analysis  proceeds  m  the  usual 
fashion.  The  1DDAM  inputs  at  each  support  for  the  first  mode  are 
obtained  from  us  weight  (1010  pounds),  its  frequency  <44  Hertz i.  and 
the  location  of  each  support  relative  to  the  point  of  impact  of  the 
torpedo.  Multiply  each  input  by  ns  factor  from  the  first  row  of  the 
participation-factor  matrix  and  sum  to  find  the  acceleration  of  the 
mode.  Scale  the  first  column  of  the  raodeshape  matrix  by  that 
acceleration  to  obtain  the  acceleration  in  the  bouncing  mode  of 
response  and  continue  to  calculate  forces,  stresses,  and  deformations 

Calculations  for  the  second  mode  will  be  similar  except  that  its 
accelerations  will  be  rotational  accelerations  m  units  of  386  radians 
per  second  squared.  Multiply  by  the  mass  moment  of  inertia  in 
pound-inches  squared  to  obtain  moment  in  pound-inches,  and 
continue  to  calculate  stresses  and  deformations  from  the  applied 
moment. 

Finally,  the  NRL  Sum  is  used  to  combine  stresses  from  the  two 
modes  with  warping  stresses  and  operating  stresses  to  estimate  peak 
stresses  at  critical  points.  The  particular  structure  of  the  example  had 
neither  warping  nor  operating  stresses  Checks  were  made  for  tensile 
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loads  m  the  bolts,  stresses  in  the  flexible  plate,  and  stresses  in  the 
welded  brackets. 


The  example  here  shows  that  the  formalism  applies  as  well  to 
rotational  responses  as  it  does  to  translations,  with  the  rotations  given 
in  radians,  stiffnesses  represented  by  moments  per  radian,  and 
masses  represented  by  mass  moments  of  inertia. 

Alternate  models  using  translations 


Figure  4  shows  four  ways  in  which  the  two  coordinates  (one 
translation  and  one  rotation)  for  the  example  could  be  replaced  by  a 
pair  of  translational  coordinates.  In  each  case  the  total  mass  of  1010 
pounds  is  represented  by  point  masses  located  so  as  to  reproduce  the 
mass  moment  of  inertia  about  an  axis  through  the  center  of  the  plate 
The  case  marked  (a)  is  of  special  interest,  since  it  shows  that 
placing  two  weights  of  150  pounds  each  directly  over  the  support 
points  not  only  gives  the  correct  mass  moment  of  inertia  but  provides 
a  total  of  278  pounds  to  match  the  effective  mass  of  the  rocking  mode. 

Each  of  the  configurations  in  Figure  4  can  be  analyzed  to 
produce  results  identical  to  the  results  obtained  using  a  rotational 
coordinate. 

A  REAL  EXAMPLE 


>ource 


Figure  5  is  adapted  from  a  sketch  in  Appendix  E  of  a  report  of 
engineering  analyses  of  submarine  equipment  done  by  a  Navy 
shipyard.1  The  figure  shows  a  4 -cubic-foot  airfiask  supported 
between  two  hull  stiffeners  by  four  beams  built  of  quarter-inch  steel 
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plate.  Calculations  from  the  analysis  are  summarized  here  witn  some 
changes  in  notation  and  revisions  in  procedures. 

Coordinates 


The  airflask  was  taken  to  have  three  degrees  of  freedom, 
represented  by  an  athwartship  translation  of  the  center  of  gravity 
f x  1 ).  rotation  about  a  vertical  axis  (xji.  and  rotation  about  a  fore-and- 
aft  axis  (X3). 

Mass  matrix 

Masses  were  represented  in  units  of  586  pounds  (slinchesi,  as 
required  by  the  computer  program  that  was  to  be  used  to  find  normal 
modes.  The  mass  for  the  translational  motion  was  117  slinches.  The 
mass  moment  of  inertia  about  the  vertical  axis  (along  the  axis  of  the 
flask)  was  taken  as  that  for  a  thin  circular  shell  with  radius  9  inches, 
or 


1.17  (9)2  -  95  lb-in-s*.  (79) 

The  shell  was  taken  as  A2  inches  long  for  its  moment  of  inertia  about 
the  fore-and-aft  axis,  or 

6(9)2  .  4(42)2 

!  17  -  [2^Z)  -  I  17i21'2  -  219  lb-m-s2  .80' 

The  mass  matrix  was  thus 

'1.17  0  O' 

M  -  ;  0  95  0  :.  '81  1 

L  0  0  2I9J 


Beams 


The  two  beams  connecting  to  Frame  90  were  angles  made  of 
quarter-inch  plate  4  inches  long,  3  inches  wide,  and  6  inches  deep 
The  cross-sectional  area  was  2  28  square  inches  and  the  area  moment 


about  the  neutral  axis  was  9  4  in4  Their  bending  stiffness  as 
cantilevers  fixed  to  the  frame  and  pinned  to  the  flask  would  be 


„  30E6(9.4) 


13  2  million  lb/in. 


Shear  stiffness  (neglecting  any  shear  factor)  would  be 


2.25(1  1E6) 


6.2  million  lb/in. 


giving  a  combined  stiffness  of 


13.2(6.2) 
13.2  *  6.2 


=  4.2  million  ib/in 


for  each  beam. 

The  two  beams  connected  to  Frame  89  were  12  inches  long,  3 
inches  wide,  and  1 1  inches  deep,  with  area  3  5  square  inches  and  an 
area  moment  of  46.4  in4  The  combined  stiffness  in  bending  and  shear 


2.413.2) 
2  4  *  3  2 


4  million  ib/in 


All  four  beams  act  against  an  athwartship  displacement  of  the 
flask,  producing  a  stiffness  of  112  million  pounds  per  inch  for 
Coordinate  1  The  beams  are  attached  7  inches  forward  and  7  inches 
aft  of  the  centerline  of  the  flask  so  that  when  the  flask  is  rotated 
about  a  vertical  axis  they  applv  a  moment  of  548  8  million  pound - 
inches  per  radian  to  Coordinate  2  Moment  arms  for  rotation  about  a 
fore-and-aft  axis  are  12  inches  above  and  below  the  center  of  the 
flask  for  a  stiffness  of  1612  8  million  pound-inches  per  radian  in 
Coordinate  3 


In  addition,  when  the  flask  is  displaced  to  starboard  (out  of  the 
paper  in  Figure  5)  the  beams  will  apply  a  moment  about  the  vertical 
axis  of 

(-4.2-4.2M7)*  (-1.4-1.4K-7)  =  -39.2  million  ib-in/in  (86) 

that  must  be  reacted  by  an  equal  and  opposite  moment  to  keep  the 
flask  from  rotating.  Similarly,  rotation  about  the  vertical  axis  will 
produce  a  net  force  per  angle  of  39.2  million  pounds  per  radian  to  port 
that  must  be  reacted  by  an  equal  force  to  starboard  to  keep  the  flask 
from  translating.  The  stiffness  matrix  needs  to  be  written  as 

f  1 1 .2E6  39.2E6  0  1 

K  =»  I39.2E6  548.8E6  0  (87) 

[  0  0  1 6 1 2.SE6J 

Support  stiffness 

A  displacement  of  Frame  90  in  the  starboard  direction  at  the 
upper  support  ( u i )  will  apply  4.2  million  pounds  of  force  per  inch  of 
displacement  to  Coordinate  1.  The  force  acts  through  a  moment  arm 
of  7  inches  to  apply  2 9.4  pound-inches  of  moment  per  inch  of 
displacement  to  Coordinate  2.  and  acts  through  a  moment  arm  of  12 
inches  to  apply  50.4  million  pound-inches  per  inch  to  Coordinate  3- 

Similar  calculations  for  the  other  three  supports  lead  to  a 
support-stiffness  matrix 

[*  4.2E6  4.2E6  1  4E6  1.4E6 

F  -  !  29.4E6  29.4E6  -9.8E6  -9.8E6  (88) 

l50.4E6  -50.4E6  16.8E6  -16.8E6. 


Normal  modes 

Data  were  entered  into  Computer  Program  RGGG1  using  the 
special  format  it  required.  The  program  returned  three  modes  with 
natural  frequencies  298.  424,  and  536  Hertz.  It  provided  a 
modeshape  matrix. 


>89) 


os 


r  1 .  o.  l.l 

1 0.185  0.  -0.066 !. 
L  0.  1.  0.  j 


that  had  been  scaled  to  make  the  largest  element  in  each  column  have 
a  value  of  unity. 

Program  RGGG1  did  not  display  the  scale  factors  that  it  used. 
They  were  recalculated  by  hand  in  a  tabulation  scheme  that  was 
equivalent  to 

r4.42  0  0  1 

S2  -  (QS)t  M  OS  -  0  219  0  ( 90 ) 

L  0  0  1 . 5  8  J 
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The  original  analysis  continued  with  rounded  calculations  to  remove 
the  squares  of  the  scale  factors  from  individual  results  as  thev 
appeared.  The  following  calculations  depart  from  that  method  to 
remove  the  scale  factors  from  the  modeshape  matrix  directly  by- 
dividing  each  column  by  the  square  root  of  its  element  of  S2.  to  obtain 

TO. 476  0.000  0.796  ] 

Q  »  i  0.088  0.000  -0.053 '.  (9!  > 

[0.000  0  068  0.000  j 

Alternate  calculation 


A  recalculation  of  the  normal  modes  was  made  using  a  1 0 -digit 
pocket  calculator  to  obtain  a  mass-normalized  modeshape  matrix 
directly  as 


■  0.483  0.000  -  0.788  ! 

Q  -  ;  -0.087  0.000  0.054  1 9 2  • 

L  0.000  0.068  0.000 j 

Differences  in  the  third  significant  figure  between  this  modeshape 
matrix  and  the  rescaled  matrix  from  RGGG1  appear  to  result  from 


w 
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rounding  of  the  scaled  modeshapes  in  Equation  89..  Frequencies  of  the 
modes  are  given  as  298.  432.  and  348  Hertz  from  the  revised  analysis. 

The  signs  are  reversed  for  rotation  about  the  vertical  axis  in 
Equation  92  because  the  data  were  entered  into  RGGG1  using  a  left- 
handed  coordinate  system.  (Direction  x  pointed  aft.  y  upward,  and  z  to 
starboard.) 

The  modeshape  matrix  from  Equation  92  has  been  used  in  the 
following  calculations. 


The  participation  factor  is 


Rl  -  P-2  QT  F 


r -0.1 56  -0.136  0.438  0.438  ; 
;  0.462  -0.462  0.154  -0.154 
0.412  0.412  0.049  0.049  j 


Sum  each  row  of  the  participation-factor  matrix  by  absolute 
value  and  square  the  sums  to  find  modal  masses  in  the  units  of  386 
pounds  used  in  the  mass  matrix.  Multiply  by  386  to  find  modal 
weights  of  550  pounds  for  Mode  1,  587  pounds  for  Mode  2.  and  328 
pounds  for  Mode  3 


One  of  the  IDDA.M  tables,  used  for  the  sake  of  an  example, 
shows  that,  at  the  point  of  impact,  a  weight  of  550  pounds  mounted  at 
298  Hertz  would  respond  with  a  peak  acceleration  of  1264  g.  A 
weight  of  587  pounds  at  432  Hertz  would  respond  to  1431  g.  and  328 
pounds  at  548  Hertz  would  respond  to  2297  g. 

One  frame  away  from  the  impact  point  the  responses  would  be 
0.45  times  as  large,  and  at  a  distance  of  24  inches  around  a  frame 
from  the  impact  point  the  responses  would  be  0.60  times  as  large 


If  a  torpedo  were  to  strike  in  way  of  the  upper  support  at  the 
after  frame  (support  uj  I.  the  IDDAM  inputs  for  the  four  supports, 
multiplied  by  the  participation  factors,  predict  peak  accelerations  for 
the  three  modes  of 


tr84  0 

o  ; 

B  -  1  0  304 

o  ;■ 
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1 5  96j 

The  scaled  raodeshape 

matrix  is 

r«o  0  1259' 

A  -  0  B  -  i  -7  0  86  I .  (95 1 

0  21  0  j 


with  accelerations  in  g  for  the  translational  coordinate  (first  row  of  the 
matrix)  and  angular  accelerations  in  multiples  of  586  radians  per 
second  per  second  for  the  rotational  coordinates  (second  and  third 
rows). 

Forces  and  moments 
Multiply 

'18  0  568  ' 

G  -  0  386  M  A  -  -269  0  3136  1 96 1 

0  1739  0  J 

to  obtain  peak  forces  (thousands  of  pounds)  applied  to  Coordinate  1 
and  peak  moments  (thousands  of  pound-inches1  applied  to 
Coordinates  2  and  3  in  each  of  the  three  modes  (The  factor  of  0  586  is 
needed  to  return  the  mass  matrix  to  convenient  units  > 


Static  analysis 

The  force  and  moment  from  Mode  1  place  a  load  of  -5  thousand 
pounds  on  each  of  the  two  after  supports,  and  a  load  of  14  thousand 
pounds  on  each  of  the  forward  supports. 

The  moment  from  Mode  2  applies  a  load  of  36  thousand  pounds 
to  each  of  the  upper  supports  and  -36  thousand  pounds  to  each  of  the 
lower  supports. 

Mode  3  puts  loads  of  254  thousand  pounds  on  each  after 
support  and  30  thousand  pounds  on  each  forward  one 

Warning 

Suppose  that  the  warping  formula  gives  a  displacement  0  200 
inch  at  the  impact  point  lui  i,  0.120  inch  at  U£,  0  090  inch  at  uj.  and 
0.054  inch  at  u«.  The  static  formula  for  the  motion  of  the  airflask  can 
be  solved  most  easily  by  writing  it  in  the  form 

x  »  K1  F  u  -  Q  P-2  QTi  F  u  =  Q  RT  u  (97  i 

that  avoids  inverting  the  stiffness  matrix  The  equilibrium  positions 
of  the  coordinates  are 

'0.116000; 

x  =  1  0.006286  i i,98  > 

L0  002875: 

The  displacement  (translation  and  two  rotations'  at  support  u\ 
is  0  194502  inch,  deflecting  the  beam  there  by  0  005498  inch  and 
producing  a  force  of  23  thousand  pounds  Similar  calculations  show 
balancing  forces  of  -25  thousand  pounds  at  U2.  -23  thousand  pounds 
at  ut.  and  *23  thousand  pounds  at  U4 


The  outermost  fibers  of  the  beams  at  Frame  90  were  4  inches 
from  the  neutral  axis  and  would  have  stress 


4  ^7  -  1  702  (99) 

psi  in  bending  per  pound  of  force  applied  at  the  end  of  the  beam. 
Shear  stress  would  be 


»  0.400 


psi  per  pound. 

The  beams  at  Frame  89  had  their  outermost  fibers  6  inches 
from  the  neutral  axis  with  stress  1  733  psi  per  pound  there  in 
bending,  and  with  shear  stress  0  286  psi  per  pound. 


The  upper  ait  beam  f u i )  has  loads,  in  thousands  of  pounds,  of  -3 
(Mode  1),  36  (Mode  2).  234  (Mode  3).  and  23  <  warping  •.  Bending 
stresses  are  -9,  62.  433.  and  39  thousand  pounds  per  square  inch.  The 
NRL  Sum  of  the  bending  stresses  is 


433  -V -9 ,2  .  1 62  ,2 


307  thousand  psi. 


i  ]  0  J } 


Shear  stresses  are  -2.  14.  102,  and  9  thousand  pounds  per 
square  inch,  with  an  NRL  Sum  of  1  1  9  thousand  psi 

Check 

The  calculation  might  as  well  stop  here  <  as  it  did  in  the  original  ' 
because  there  is  no  chance  that  the  support  at  the  impacted  frame 
could  withstand  the  bending  stress  of  507  thousand  pounds  per 
square  inch  or  the  shear  stress  of  1  19  thousand  pounds  per  square 
inch  The  bending  stress  is  more  than  ten  times  the  yield  stress  for  a 
medium  steel  and  the  shear  stress  is  more  than  four  times  yield  in 
shear.  A  redesign  of  the  supports  for  the  airflask  is  obviously  needed 
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APPENDIX  A 

A  SHORT  INTRODUCTION  TO  MATRIX  ALGEBRA 

MATRICES 

A  computer  statement  similar  to 

DIMENSION  M(4.5>  »A-H 

will  reserve  20  memory  spaces  for  numbers  that  can  thereafter  be 
referenced  as  M(J,L),  with  J  an  integer  in  the  range  from  1  to  4  and  L 
an  integer  in  the  range  from  1  to  5  The  numbers  are  usually  thought 
of  as  arranged  in  a  rectangular  array  with  four  rows  and  five  columns, 
designated  simply  by  the  letter  M. 


"Ml  1,1)  M(  1 ,2 
MI 2,1)  M( 2.2 


|  Ml 3,: 
LM<  4. : 


Mil, 2)  M(  1 .3)  MU .4)  M ( 1.5) 
Ml 2.2)  MI2.3')  MI 2.4)  MI2.5) 
M(3.2)  M( 3.3 )  MI 3.4)  MI3.5) 
MI 4.2)  MI  4,3)  MI 4,4 1  M14.5) 


SOME  DEFINITIONS 

The  matrix  with  dimension  MI 4.3 1  is  called  a  rectangular  matrix 
because  it  has  a  different  number  of  rows  than  columns.  A  matrix 
with  dimension  MI 4,4)  would  be  called  a  square  maim  The 
elements  MU.J)  are  called  the  diagonal  elements .  they  would  be  the 
elements  M'  1,1  f  MI2.2),  MI 3.3  h  and  M'  4.4)  for  the  matrix  here 

The  sum  of  the  diagonal  elements  is  called  the  trace  of  the 
matrix;  the  trace  of  the  matrix  here  would  be  MI  1,1 )  *  M',2.2)  -  M1 3,3  1 
-  MI  4.4  i  A  symmetric  maim  is  one  in  which  everv  MiJ.L)  is  equal  to 
the  corresponding  M<L.J)  -  that  is  Mi  1.2'  =  Ml 2. If  Mt  1.3)  =  M '3.1  1 
and  so  on  for  all  the  pairs 

SOME  OPERATIONS 

If  matrix  0  has  elements  QlJ.Lf  the  matrix  with  element?  L.j ' 
is  called  the  transpose  of  Q  designated  by  Transposition 

corresponds  to  interchanging  rows  and  columns  in  the  arrav  If  the 


m 


transposed  matrix  is  equal  to  the  original  matrix  tQT  -  Qj,  then  the 
matrix  must  have  been  both  square  and  symmetric. 

Matrices  of  the  same  dimension  are  added  or  subtracted  by 
adding  or  subtracting  the  corresponding  elements.  Matrices  of 
different  dimensions  cannot  be  added  to  or  subtracted  from  one 
another. 

MATRIX  MULTIPLICATION 

The  product  of  two  matrices  is  written  as  if  it  were  ordinary 
multiplication.  That  is. 

Z-  MO  (A -3) 

is  the  product  of  M  and  Q.  The  operation  is  not  an  element-by¬ 
element  multiplication,  however,  but  involves  multiplying  each 
element  in  a  row  of  the  first  matrix  by  the  corresponding  element  in  a 
column  of  the  second  matrix  and  summing  the  products  Each  element 
of  Z  is  generated  as 

Z(J,N)  -  M(J.l)  Q(  1.N)  -  M(  1.2)  Q(2,N)  *  M(  1.3)  Q«3.N) 

*  M(J.4)Q(4.N)  *  ....  (A -4) 

The  equation  can  be  abbreviated  as 

Z(J.N)  =  M(J,L)  QlL.Ni.  i  A  -  5  ‘ 

where  it  is  understood  (the  summation  convention  i  that  the  product 
is  to  be  summed  over  ail  possible  values  of  the  repeated  index  L. 

The  product  matrix  Z  will  have  as  many  rows  as  there  were 
rows  in  M,  and  as  many  columns  as  there  were  columns  in  Q 
Moreover.  M  must  have  exactly  a?  many  columns  as  Q  has  rows,  so 
that  the  index  L  can  be  assigned  to  both  and  run  through  a  single 
range  of  values 

Because  of  the  rule  for  multiplication,  it  is  not  surprising  that 
the  matrix  product  QM  is  a  different  matrix  from  the  product  MQ 
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MORE  DEFINITIONS 


A  diagonal  matrix  is  one  in  which  all  of  the  elements  other  than 
the  diagonal  elements  have  values  of  zero.  A  untt  matrix  (usually 
designated  by  I)  is  a  diagonal  matrix  in  which  each  diagonal  element 
has  a  value  of  unity. 

The  inverse  of  a  matrix  K  is  the  matrix  K'1  that  makes 

K1  K  =  I  (A -6) 

be  a  unit  matrix.  Not  all  matrices  have  inverses.  The  definition  of  an 
inverse  is  reciprocal;  K  is  also  the  inverse  of  K*1  so  that 

K  K'1  =  I.  ( A -7 ) 

If  the  transpose  of  a  matrix  is  also  its  inverse,  the  matrix  is  said 
to  be  unitary  That  is. 

UTU  *  I  (A -Si 

defines  U  as  a  unitary  matrix.  A  simple  example  is 
r  cos(  H )  -sint  H  T; 

U  "  -Sin(H)  cos'  H )  J‘  l  A -9  i 

which  is  unitary  for  any  angle  H. 

VECTORS 

A  matrix  with  only  one  column  is  called  a  column  vector  If  a 
matrix  X  has  dimension  X'  5, 1 ).  for  example,  it  is  a  column  vector  with 
five  elements  It  is  usual  to  suppress  the  unit  uesignation  and  show 
its  dimension  as  X( 5 )  Column  vectors  are  designated  by  lower-case 
boldface  letters,  such  as  i.  in  the  present  report. 

A  matrix  with  onlv  one  row  is  called  a  row  vector  h  is 
dimensioned  X*  1.5  c  for  example,  and  is  designated  here  using  the 
symbol  for  transpose.  xT  Again  the  unit  designation  is  usuallv 
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suppressed  to  show  a  row  vector  as  dimension  X(5).  although  this  has 
the  defect  of  not  distinguishing  it  from  a  column  vector. 

All  of  the  rules  for  adding,  subtracting,  and  multiplying 
matrices  apply  to  vectors  provided  the  unit  dimension  is  understood. 

SOME  HELPFUL  RELATIONS 

The  following  rules  can  be  verified  by  working  out  some  simple 
examples.  First,  the  transpose  of  a  matrix  product  is  obtained  from 
the  transposes  of  the  factors  taken  in  the  reverse  order,  or 

(MQ)t  =  Qt  MT  1 A  - 1 0 1 

Similarly,  the  inverse  of  a  product  involves  inverting  and 
interchanging  the  factors,  as 

(MQ)'1  =■  Q-l  M-l  !  A  - 1 1  i 

If  a  matrix  is  symmetric,  its  inverse  is  also  symmetric. 
However,  the  product  of  two  symmetric  matrices,  such  as  M  1  and  K. 
is  not  a  symmetric  matrix,  since 

(M-l  K)t  -  KT  M-T  -  KM'1  I A  - 1 2 1 

is  different  from  the  original  product  M  1  K 

Both  the  mass  and  stiffness  matrices  M  and  K  for  a  structure 
are  necessarily  symmetric  If  a  stiffness  matrix  were  not  symmetric, 
for  example,  the  structure  could  be  run  as  a  perpetual-motion 
machine  by  moving  it  around  a  cycle  in  which  the  smaller  off-diagonal 
element  was  worked  during  the  compression  phase  and  the  larger 
offdiagonal  element  was  used  to  provide  power  during  the  return 
stroke 

Matrix  elements  for  a  structure  can  often  by  generated  or 
checked  by  using  the  useful  relations  that  K  x/2  is  the  potential  'or 
strain)  energy  and  M  v/2  is  the  kinetic  energy,  where  x  is  a 
column  vector  of  displacements  and  v  a  column  vector  of  velocities 


USEFULNESS  OF  MATRIX  ANALYSIS 


The  rules  for  addition,  subtraction,  and  multiplication  of 
matrices  were  chosen  specifically  so  that  matrix  equations  very 
similar  to  scalar  equations,  such  as 

Ma  +  Ki=Fu,  l  A- 13) 

would  expand  into  precisely  the  set  of  coupled  differential  equations 
needed  to  represent  a  structure  with  multiple  coordinates  The  great 
compaction  of  the  notation  and  the  resultant  improvement  in 
comprehension  are  exactly  what  is  needed  in  analyzing  complicated 
structures. 


APPENDIX  B 

A  COMPUTER  PROGRAM  FOR  RIDDAM 

PROGRAM 

1000  REH  THIS  PROGRftn  IS  WRITTEN  IN  CORE  BRSIC 

1100  REN  I  DENT  I FV  THE  PROGRfiN 
1110  PRINT  "877F7  RIOOfln" 

1120  PRINT  "  ACCEPTS  fi  0 i RGONRL  NRSS  flRTR I X  IN" 

1130  PRINT  "  POUNDS  OR  POUND  INCHES  SOURRED." 

1H0  PRINT  “  STIFFNESSES  IN  POUNDS  PER  INCH  OR" 

1150  PRINT  “  POUND- INCHES  PER  RRDIfiN," 

1160  PRINT  “R  L  SORT,  NRURL  RESEARCH  LAB,  UfiSH I NGTOM ,  DC 


1170  PRINT  " JULV  15,  1987." 

1180  PRINT 

1200  REN  D I  HENS  I  ON  THE  URRiflBLES 

•:,r’  :,M  M/,c,:o),  ko: 

1220  DIN  0(10,10),  P(10, 10),  R( 10, 10) 

1230  DIN  FI ( 1 0) ,  U(10),  0(10),  F2( 1 0) , B( 1 0) 

12- 10  DIN  fl(10, 10), GOO,  10) 

1300  REN  CASE  NUNBER  OR  TEST  CASE 

1310  PRINT  "ENTER  fl  NUNBER  FOR  VOUR  CRSE,  OR" 

1320  PRINT  "ENTER  0  (ZERO)  TO  RUN  fi  TEST  CfiSE." 

1330  INPUT  T 

13- 10  IF  T >0  THEN  MIQ 

1350  PRINT  "CASE  0  IS  fi  IC'OO-POUND  BAR  SUPPORTED" 

1360  PRINT  “  INDEPENDENTLY  AT  EACH  END.  IT  IS  NODELED  AS" 
1370  PRINT  "N1  EQUAL  NASSES  UlTH  Ml -1  EQUAL  SPRINGS." 

1380  PRINT 


MOO  REn  ASSIGN  DEGREES  OF  FREEDOH 

HI 0  PRINT  "ENTER  DEGREES  OF  FREEDOH  (1  TO  10) 

1 420  INPUT  N1 

1430  IF  HI < 1  THEN  1410 

1440  IF  N1 >10  THEN  1410 

1500  REfl  CLERR  SONS  HRTRICES 

1510  FOR  J-1  TO  N1 

1520  FOR  L-1  TO  HI 

1530  LET  n(J,L)«0 

1540  LET  K(J,L)=0 

1550  LET  Q(J,L)-0 

1560  NEXT  L 

1570  NEXT  J 

1580  IF  T >0  THEN  1010 

1600  REH  TEST  CASE,  HASS  AND  STIFFNESS 
1610  LET  N2=2 
1620  FOR  J-1  TO  N1 

•  WOVJ  Lwi  IIWjJ/2!  JvJU/  v,  i  *  »  J 

1640  LET  K(J(J)-2*20000*(NM1) 

1650  NEXT  J 

1660  FOR  J-1  TO  Ml -1 

1670  LET  K( J, J+1 )«-20000*(N1 *1 ) 

1680  LET  K(  J+1 , J)=K( J, J+1 ) 

1690  NEXT  J 

1700  REH  TEST  CASE,  SUPPORT  STIFFNESS 

1710  FOR  J-1  TO  Nl 

1720  LET  F( J, I )=0 

1730  LET  F(J,2)=0 

1740  NEXT  J 

1750  LET  F( 1 , 1 )*20000*(N1 +1 ) 

1760  LET  F ( N 1 ,2)*20000*(N1  *1 ) 

1770  GOTO  2410 


1800  REn  ASSIGN  NUflBER  OF  SUPPORTS 

1810  PRINT  "ENTER  NUflBER  OF  SUPPORTS  (1  TO  10)" 

1820  INPUT  N2 

1830  IF  N2< 1  THEN  1810 

1840  IF  H2 > 1 0  THEN  1810 

2000  REN  ENTER  THE  HATRIX  ELEHENTS 

2100  REH  DIRGONRL  HASS  DfiTR I X  n 
2110  FOR  J-1  TO  N1 
2120  PRINT  "ENTER  nfiSS" ; J 
2130  INPUT  n(J,J) 

2140  IF  n(J,J)>0  THEN  2170 

2150  PRINT  "nflSSES  HUST  5E  POS I T :  'JE  .  "PV  AGAIN. 
2160  GOTO  2120 
2170  NEXT  J 

2200  REn  SVnnETRlC  STIFFNESS  nfiTP,,x  c 
2210  FOR  J-1  TO  N1 
2220  FOR  L=J  TO  N1 

2230  PRINT  "STIFFNESS  HATR • X  FCIT.  ;  "SOL jHN"  l 
2240  INPUT  < ( J , L ) 

2250  LET  < ( L , J )  =< 1 J , L ) 

2260  NEXT  L 
2270  NEXT  J 

2300  REn  SUPPCRT-ST . FFNt:  ;  •  - 

2310  FOR  j*l  TO  n ' 

2320  FOR  L* 1  tq  n: 

2330  PRINT  -cy . OOORC  '•"=*-■  ""  ;  pc-ct 

2310  INP!JT  F  :  L 
2350  NEXT  . 

2360  NEXT  j 
2370  PRINT 


2400  REM  LABEL  THE  CASE 
2410  PRINT  “THIS  IS  CASE“;T 
2420  PRINT 

2500  REH  ECHO  THE  HASS  MATRIX 
2510  PRINT  “MASS  MATRIX  M  =" 

2520  FOR  J-1  TO  N1 
2530  FOR  L-l  TO  N1 
2540  PRINT  H(J,L), 

2550  NEXT  L 
2560  PRINT 
2570  NEXT  J 
2580  PRINT 

2600  REM  ECHO  THE  STIFFNESS  MATRIX 
2610  PRINT  "STIFFNESS  MATRIX  <  •  “ 
2620  FOR  J-1  TO  N1 
2630  FOR  L-l  TO  Nl 
2640  PRINT  K( J , L ) , 

2650  NEXT  L 
2660  PRINT 
2670  NEXT  J 
2680  PRINT 


2?0G  REn  ECHO  THE  SUPPORT-STIFFNESS  MATRIX 
27’  0  PRINT  “SUPPORT-STIFFNESS  MATRIX  F  «“ 

2 7 2  'j  FOR  j-1  TO  Ml 
2730  FOR  L-l  TO  N2 
2 " 4 n  f ( j , l ) , 

2"5C 

'  ->.Q  DC 
'770  *CyT 
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3000  RER  CHANGE  OR  CORRECT  RATRIX  ELERENTS 
3100  REH  CHECK  THE  flfiTRICES 

3110  PRINT  "PLEASE  CHECK  THE  RATRICES.  IF  THEV  ARE  GOOD 
3120  PRINT  "ENTER  0.  IF  VOU  URNT  TO  CHANGE  ANV,  ENTER  1. 
3130  INPUT  J 
3140  IF  J=0  THEN  4110 

3200  REN  CHANGE  THE  NASS  RATRIX  n 

3210  PRINT  “RASS  TO  BE  CHANGED  OR  0  IF  NO -CHANGE" 

3220  INPUT  J 

3230  IF  J<1  THEN  3310 

3240  PRINT  f1(  J,  J) ; "  IS  OLD  nfiSS";J;",  ENTER  MEL!  UALUE ; " 
3250  INPUT  R( J, J) 

3260  IF  R(J,J)>Q  THEN  3210 

3270  PRINT  “RASSES  RUST  BE  POSIT IUE  NUMBERS ! “ 

3230  GOTO  3240 

3300  RER  CHANGE  THE  STIFFNESS  RATRIX  K 

3310  PRINT  "ENTER  ROU,  COLURN  OF  ELERENT  TO  BE" 

3320  PRINT  "CHANGED  IN  K.  ENTER  0,0  IF  NO  CHANGE." 

3330  INPUT  J, L 
3340  IF  J=0  THEN  3410 

3350  PRINT  K(J,L); "IS  OLD” ; J ; L ; " .  ENTER  NEU  UALUE:" 

3360  INPUT  K( J, L) 

3370  LET  K(L, J)=K( J, L) 

3380  GOTO  3310 


1 
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3400  REn  CHANGE  THE  SUPPORT-STIFFNESS  flRTR I X  F 
3410  PRINT  "ENTER  ROU,  COLUHN  OF  ELEHENT  TO  BE" 

3420  PRINT  "CHANGED  IN  F.  ENTER  0,0  IF  NO  CHRNGE. " 

3430  INPUT  J,L 

3440  IF  J=0  THEN  2410 

3450  PRINT  F(J,L); “ IS  OLOH;J;L;“.  ENTER  NEU  URLUE: “ 
3460  INPUT  F(J,L) 

3470  GOTO  3410 

4000  REN  START  NORNfiL-nOOE  RNRLVSIS 
4100  REfl  NORMAL  I ZE  ON  THE  MRSS  MRTR IX,  QT*n*q  =  I 
4110  FOR  J-l  TO  N1 
4120  LET  q(J>J)«1/SQR(H(JJJ)) 

4130  NEXT  J 

4200  REn  SVnnETRlC  OVNflniC  nflTRIX,  QT*K*q 
4210  FOR  J-1  TO  N1 
4220  FOR  L-1  TO  HI 

4230  LET  P(J, L)=Q( J, J)*K(J,L)*0(L, L) 

4240  NEXT  L 
4250  NEXT  J 

4300  REn  SHOU  THE  OVNRM I C  MATRIX 
4310  PRINT  "SVnnETRlC  OVNflniC  HftTR ! X  QT*K*Q  =“ 

4320  FOR  J*1  TO  HI 
4330  FOR  L«1  TO  HI 
4340  PRINT  P ( J , L ) , 

4350  NEXT  L 
4360  PRINT 
4370  NEXT  J 
4380  PRINT 


1  ^  ^ 


5000  REn  D ! flGONRL I ZE  THE  DVNfitllC  (1RTR I X 

5100  REH  SET  COUNTER  AND  RUN  THROUGH  OFFDIRGONRL  ELEHENTS 
5110  PRINT  "CALCULATING. ..." 

5120  LET  E=0 
5130  FOR  J-1  TO  N1-1 
5140  FOR  L=J+1  TO  N1 

5200  REN  SORT  COLUHNS  8V  INCREASING  FREQUENCV 
5210  IF  P( J, J)<=P(L,L)  THEN  5110 
5220  FOR  J1-1  TO  N1 
5230  LET  Z»Q(J1,J) 

5240  LET  0( J1 , J)*Q(J1 , L) 

5250  LET  0(J1.L)-Z 
5260  LET  Z=P(J1,J) 

5270  LET  P ( J 1 , J)=P(J1 fL) 

5280  LET  P(J1,L)«Z 
5290  NEXT  J1 

5300  REN  RNO  INTERCHANGE  ROUS  OF  P  fib  UELL 
5310  FOR  JI-1  TO  N! 

5320  LET  2*P( J, J1 ) 

5330  LET  P(J,J1)*P(L,J1) 

5340  LET  P(L,J1)-Z 
5350  NEXT  J1 

5400  REn  SKIP  SHALL  UALUES 

5410  IF  A6SfP(J,L)+P(L,J))< .Cl  TnEN  5920 

5500  REH  DON ‘ T  DIUIDE  BV  ZERO 
5510  IF  P(J, J ) < P ( L , L )  THEN  5610 
5520  LET  C'=SOR(  .5) 

5530  LET  S-C 
5540  GOTO  5710 
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5600  REfl  JACOB  I  ROTATION 

5610  LET  H».5*ATN((P(JjL)+P(LjJ))/(P(J)J)-P(L,L))) 
5620  LET  C*COS(H) 

5630  LET  S*S I N(H) 

5700  REN  UPDATE  THE  COLUNNS  OF  THE  HATRICES 

5710  FOR  J1-1  TO  HI 

5720  LET  Z*C*Q(J1,J)*S*Q(J1,L> 

5730  LET  Q(J1 , L)=C*Q(J1 ,L)-S*Q(J1 , J) 

5740  LET  Q( J1 , J)=2 

5750  LET  Z=C*P( J1 ,  J)*S*P(J1 ,L) 

5760  LET  PC J1 ,L)-C*P(J1 ,L)-S*P(J1 ,J) 

5770  LET  P(J1,J)-2 
5780  NEXT  J1 


5800  RED  AND  ALSO  THE  ROUS  OF  P 

5810  FOR  J1-1  TO  HI 

5820  LET  Z3C*P( J, J1 )+S*P(L, J1 ) 

5830  LET  P(L, J1 )=C*P(L, J1 )-S*P(J, J1 ) 
5810  LET  P(J,J1)-2 
5850  NEXT  J1 


5900  REN  STEP  THE  COUNTER  ANO  END  THE  LOOP 
5910  LET  E-E-1 
5920  NEXT  L 
5930  NEXT  J 


6000  REN  SHOU  AND  TELL  THE  NORHAL  NODES 


6100  REN  IS  P  0 1  AGONAL  VET7 

6110  PRINT  "ZAPPED" ;E, "OF  THE  OFFO 1  AGONAL  ELE"E*T: 
6120  IF  E  >0  THEN  5120 


6200  REH  SHOU  nOOESHRPES 

6210  PRINT  "MODESHflPE  flRTR I X  Q  =“ 

6220  FOR  J=1  TO  HI 
6230  FOR  L-l  TO  N1 
6240  PRINT  O(J.L), 

6250  NEXT  L 
6260  PRINT 
6270  NEXT  J 
6280  PRINT 

6300  REN  SHOU  FREQUENCV-SQURREO  flRTR  I X 

621 Q  PRINT  “FREQUENCV-SQURPED  MATRIX  QT*K*Q 

6320  FOR  J»1  TO  N1 

6330  FOR  L*J  TO  N1 

6340  PRINT  P(J,L), 

6350  NEXT  L 
6360  PRINT 
6370  NEXT  J 
6380  PRINT 

6100  REH  CHECK  FOR  LOOSE  PRRTS 

6110  FOR  J*1  TO  N1 

6120  IF  P(J,J)<>0  THEN  6110 

6130  LET  P(J,J)--.01 

6110  NEXT  J 

6500  REH  PORT  I C  I  PfiT I ON-FRCTOn  MATRIX  R 


6510  FOR 


TO  N 1 


6520  FOR  L-i  TO  N2 
6530  LET  R ( J , L ) =  0 
6510  FOR  J 1  *  1  TO  N! 

6550  LET  ?.( J,L)»R(J.L)*Q( Ji . J 1  *c 1  J 1 . 
6560  NEXT  ji 

6570  LET  Rf  j  L  >=R(  J  ,  L  )  /P (  J ,  j  • 

6560  NEXT  L 
6590  NEXT  „ 


6600  REH  SHOU  THE  PRRT I C I PAT I  OH  FRCTOR 
6610  PRINT  "PARTICIPATION  FACTOR  R  =" 

6620  FOR  J-1  TO  N1 
6630  FOR  L=1  TO  N2 
6610  PRINT  R ( J , L ) , 

6650  NEXT  L 
6660  PRINT 
6670  NEXT  J 
6680  PRINT 

6700  REH  CHECK  FOR  NEGRTIUE  E I GENUALUES 

6710  FOR  J-1  TO  N1 

6720  IF  P( J, J) >0  THEN  6770 

6730  IF  P( J , J) >- .011  THEN  6760 

6710  PRINT  "STRUCTURE  COLLAPSES  IN  nODE“;J 

6750  PRINT  "THERE  IS  fl  fll STAKE  IN  VOUR  STIFFNESS  flRTR I X 

6760  LET  P(J,J)=0 

6770  NEXT  J 

6800  REN  CALCULATE  FREQUENCIES 
6810  FOR  J-1  TO  N1 

6820  LET  FI(J)*SQR(386*P(J,J))/6.283185 
6830  NEXT  J 

6900  REN  CALCULATE  NODAL  UEIGHTS 

6910  FOR  J-1  TO  N1 

6920  LET  U(J)-0 

6930  FOR  L-1  TO  N2 

6910  LET  U( J)=U( J)+fiBS(R( J, L) ) 

6950  NEXT  L 

6960  LET  U(J)-U(J)*U(J) 

6970  NEXT  J 
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7000  REn  FREQUENCIES,  UEIGHTS,  RND  INPUTS 


v\ 
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7100  REH  REAL  OR  FAKE  INPUTS 
7110  IF  T=0  THEN  7110 

7120  PRINT  "ENTER  0  TO  GET  FAKE  INPUTS  (EXRNPLE  GNLV) , 
7130  PRINT  "OR  ENTER  1  TO  SUPPLV  INPUTS  FROtl  TRBLES. “ 
7110  INPUT  J 
7150  IF  J-0  THEN  7110 

7200  REH  ENTER  INPUTS 
7210  FOR  J-1  TO  N1 

7220  PRINT  "INPUT  FOR" ; U( J) ; "POUNOS  ANO" ; FI (J); "HERTZ; 
7230  INPUT  D( J) 

7210  NEXT  J 

7300  REN  ENTER  FRRNE  FACTORS 
7310  FOR  J-1  TO  N2 

7320  PRINT  "FRANE  FACTOR  FOR  SUPPORT" ;J 
7330  INPUT  F2(J) 

7310  NEXT  J 
7350  GOTO  7610 

7100  REN  INPUTS  (FAKE)  FOR  TEST  CRSE 
7110  FOR  J-1  TO  HI 

7120  LET  D( J)3* . 3*P( J, J)/$QR(  1  *.000001 *U(J)*P(J,  J) ) 
7130  NEXT  J 

7500  REN  FRflHE  FACTORS  (FAKE)  FOR  TEST  CASE 

7510  LET  F2(1 )-1 

7520  FOR  J=2  TO  N2 

7530  LET  F2(J)-.5*F2(J-1) 

7510  NEXT  J 
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7600  REH  SHOU  THE  HOOES 

7610  PRINT  "nOOEVFREQUENCY“,“UEIGHTVINPUT" 
7620  PRINT  , "(HERTZ)", "(POUNDS)”, "(G)" 

7630  FOR  J-1  TO  NI 

7640  PRINT  J,  F1(J)jU(J)i0(J) 

7650  NEXT  J 
7660  PRINT 

7700  REN  SHOU  THE  FRfiflE  FACTORS 

7710  PRINT  "FRAflE  FACTORS  (RATIO  TO  i HPRCT  POINT) 

7720  FOR  J-1  TO  N2 

7730  PRINT  F2(J) , 

7740  NEXT  J 
7750  PRINT 
7760  PRINT 

8000  REN  PEAK  RESPONSES 

8100  REM  CALCULATE  ACCELERATIONS  OF  nODES 

8110  FOR  J-1  TO  NI 

8120  LET  8( J)=Q 

8130  FOR  L»1  TO  N2 

8140  LET  B(J)»B(J)*R<J,L)*F2(L)*D(J} 

8150  NEXT  L 
8160  NEXT  J 


8200  REN  CONUERT  TO  ACCELERATIONS  OF  HASSES 

8210  FOR  J-1  TO  NI 

8220  FOR  L-1  TO  NI 

8230  LET  A( J, L)=Q( J,  L)*B(L) 

8240  NEXT  L 
8250  NEXT  J 


8300  REH  FI  NO  PEAK  FORCES  IN  EACH  NODE 

8310  FOR  J*l  TO  HI 

8320  FOR  L=1  TO  NT 

8330  LET  G(J,L)=f1(  J,  J)*R(J,L) 

8310  NEXT  L 
8350  NEXT  J 

8100  REN  SHOU  THE  ACCELERATIONS 

8110  PRINT  "PEAK  ACCELERATIONS  (6)  IN  EACH  nOOE" 

8120  FOR  J-1  TO  N1 

8130  FOR  L-1  TO  N1 

8110  PRINT  A( J,L) , 

8150  NEXT  L 
8160  PRINT 
8170  NEXT  J 
8180  PRINT 

8500  REH  SHOU  THE  FORCES 

8510  PRINT  "PEAK  FORCES  (POUNDS)  IN  EACH  NODE" 

8520  FOR  J=1  TO  N1 
8530  FOR  L=1  TO  N1 
8510  PRINT  G( J , L) , 

8550  NEXT  L 
8560  PRINT 
8570  NEXT  J 
8580  PRINT 

9000  REH  UINOUP 

9100  REn  UMRT  3HPILL  UE  00  HEMT~ 

9110  PRINT  “THIS  IS  ALL  THIS  PROGRAn  DOES.  CHOOSE 
9120  PRINT  "1  -  NEU  IflPACT  PO I  NT- 
9130  PRINT  "2  -  NEU  INPUT  TABLES" 

9110  PRINT  "3  -  CHANGE  THE  SUPPORTS" 

9150  PRINT  "1  -  CHANGE  THE  STIFFNESS  nflTRiX" 

9160  PRINT  “5  -  CHANGE  THE  HASS  HATRIX" 

9170  PRINT  "6  -  NEU  STRUCTURE" 

9180  PRINT  "9  -  END  THIS  RUN" 


9200  REfl  CHOI 
9210  INPUT  J 
9220  IF  J=1  1 
9230  IF  J=2  1 
9240  IF  J=3  1 
9250  IF  J=4  1 
9260  IF  J=5  1 
9270  IF  J-6  1 
9280  IF  J=9  1 
9290  GOTO  911 


THEN  7310 
THEN  7210 
THEN  3410 
THEN  3310 
THEN  3210 
THEN  1310 
THEN  9910 
110 


9900  RED  END 

9910  PRINT  “THANK  YOU  FOR  RUNNING  THIS  PROGRRh." 

9920  PRINT  "FINISHED" 

9930  END 

COMMENTS 

Documentation 

The  program  is  self -documented  by  remarks.  !;  uses  the 
same  procedures  and  most  of  the  same  svmbols  a?  are 
described  m  the  text  of  this  report 

Core  BASIC 

The  program  is  written  in  primitive  DAS  IT  .uncage 
without  any  of  the  enhancements  that  are  available  m  var.>.us 
dialects  of  BASIC.  As  such,  it  is  long  and  slow  but  should  run  as 
is  on  any  computer  that  has  a  BASIC  compiler  or  ;n'.erp^e'.er 
The  program  is  also  wasteful  of  memory  since  f;  r  c.ar-i\  ;t 
does  not  use  anv  shortcuts  and  defines  new.  variables  or  earn 
step  in  the  calculation 
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SAMPLE  OUTPUT 

Test  case.  Case  0 

The  program  is  loaded  with  an  example  consisting  of  a 
unilorm  bar  with  weight  1000  pounds  and  axial  stiffness 
20,000  pounds  per  inch.  The  bar  can  be  represented  by  from 
one  to  ten  lumped  masses  evenly  spaced  along  a  massless 
spring  A  simple  analytical  solution  is  available  for  the  fixed- 
base  frequencies  and  modeshapes  of  the  string  of  equal 
masses  - 

Following  is  part  of  the  output  from  Case  0  The  bar  was 
modeled  using  four  stations,  each  representing  the  motion  of  a 
200-pound  mass  lump  The  remaining  mass  was  lumped  with 
the  supports  at  each  end  of  the  bar 

The  test  case  uses  a  simple  formula  to  generate  IDDAM- 
hke  inputs  The  inputs  are  for  sake  of  the  example  only  and  are 
not  intended  to  applv  to  anv  particular  situation. 
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SUPPORT-STIFFNESS  MATRIX  F  = 

1 OOOOO  0 

0  0 

0  0 

0  100000 

svnnETRic  dynamic  hrtr i x  qt*k*q  = 

1000  -500  0 

-500  1000  -500 

0  -500  1000 

0  0  -500 

CALCULATING. . . . 

ZAPPED  6  OF  THE  OFFO I  AGONAL  ELEMENTS. 

ZAPPED  6  OF  THE  OFFD I  AGONAL  ELEMENTS . 

ZAPPED  5  OF  THE  OFFD  I  AGONAL  ELEflENTS . 

ZAPPED  0  OF  THE  OFFO  I  AGONAL  ELEflENTS. 


0 

0 

-500 

1000 


flODESHAPE  flRTR  I X  Q  = 
2.623614E-02  -4.25323E-02 
4.253229E-02  -2.626693E-02 


4.25323E-02 

2.628698E-02 


2 . 62861  IE-02 
4 . 253229E-02 


4.253267E-G2 

-2.628675E-02 

-2.628635E-02 

4.253241E-02 


FREOUENCV-SQUARED  MATRIX  QT*K*Q  * 

190.983  4 . 936953E-03  2.26569E-05 

4.9I2365E-03  690.983  -5.086626E-06 

3 . 227984E-05  1.204654E-06  1309.017 

-7 . 09665E-06  2.624094E-05  -2 .425924E-03 

PARTICIPATION  FACTOR  R  = 

13.7636  13.76404 


-2.628634E-02 
4 . 25324E-02 
-4 , 253267E-02 
2 . 628677E-02 


5 . 789629E-06 
•2.625585E-05 
-2 . 479627E-03 
1809.017 


-6.155405 
3.249208 
-1 .453074 


13.76404 
6 . 1 55331 
3.249187 
1.453097 


